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ABSTRACT 


A  near  miss  underwater  explosion  can  significantly  damage 
improperly  shock  hardened  combat  systems  equipment  and  render 
the  ship  uncible  to  “fight  hurt*'.  MIL-S-901D  currently  requires 
shock  qualifying  mediximweight  equipment  to  a  "generic"  shock 
excitation  on  the  Navy's  Mediumweight  Shock  Machine  (MWSM) . 
This  shock  excitation  is  severe,  but  not  always  characteristic 
of  the  actual  ship  structure  response  to  an  underwater 
explosion.  This  study  proposes  a  design  modification  which 
will  allow  using  a  multi-DOF  equipment  mounting  fixture  on  the 
MWSM  which  can  be  "tuned"  to  simulate  shipboard  shock 
characteristics  determined  from  modal  testing  or  previous  ship 
shock  trial  data.  Equipment  qualified  in  this  manner  could  be 
highly  relied  on  to  survive  in  battle. 
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1. 


INTRODUCTION 


In  today's  world  of  modern  warfare  the  U.S.  Navy  must  rely 
on  state  of  the  art  combat  systems  equipment  to  maintain  its 
fighting  advantage  over  possible  advisaries .  However,  the 
high-tech  nature  of  this  equipment  can  lead  to  increased 
vulnerability  to  mechanical  shock  induced  failure  if  not 
adequately  packaged  to  withstand  the  severe  excitations 
expected  from  conventional  and  nuclear  underwater  weapons . 
When  exploded  in  proximity  of  the  ship,  these  weapons  produce 
an  intense  pressure  wave  v;hich  is  applied  over  the  entire 
underwater  portion  of  the  ship's  hull.  Although  the  hull  is 
likely  to  remain  intact,  the  violent  complex  shock  waves  that 
propagate  throughout  the  ship  can  sufficiently  damage 
essential  equipment  to  render  the  ship  unable  to  "fight  hurt" 
and  significantly  inpair  its  mission  integrity. 

The  Navy  currently  uses  Military  Specifications 
(MIL-S-901D) ,  “Shock  Tests,  High  Impact;  Shipboard  Machinery, 
Equipment  6uid  Systems,  Requirements  For"  to  detail  specific 
shock  qualification  reejuirements  for  shipboard  machinery, 
equipment,  systems  and  structures  which  are  required  to  resist 
the  effects  of  mechanical  shock.  These  requireme.nts  establish 
the  general  shock  test  criteria  and  provide  the  contracting 
activity  a  basis  for  selecting  the  appropriate  testing  device 
based  on  the  weight  category  of  the  ecfuipment.  The  weight 
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rjategories  defined  are;  lightweighc  for  an  attached  weight  up 
to  550  lb;  mediumweight  for  an  attached  weight  up  to  7400  lb; 
and  heavyweight  for  a  total  weight  up  to  60,000  lb.  The 
lightweight  and  mediumweight  machines  are  similar  in  that  the 
high  impact,  shock  is  delivered  to  the  attached  equipment  by 
use  of  a  hammer  and  amvil  assembly.  Hea'/yweight  category 
equipment  is  installed  onboard  a  floating  platform  barge  and 
subjected  to  shock  from  an  underwater  explosive  of  known 
charge  and  standoff  geometry.  Of  the  three  weight  categories 
listed,  this  research  focuses  on  the  use  of  the  Medium  Weight 
Shock  Machine  (MWSM)  shown  in  Figure  1 . 


Figure  1.  The  havy  High-Iitpact  Shock  Machine  for 
Mediumweight  Equipment  {MWSM).  Courtesy  of  Clements (1972 ) . 
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The  MWSM  delivers  a  vertical  high  impact  mechanical  shock 
to  the  anvil  table  by  means  of  a  3000  lb  swinging  hammer.  The 
impact  between  the  hamm.er  auid  anvil  is  highly  elastic  and  the 
energy  is  controlled  by  adjusting  the  height  of  the  hammer 
above  the  anvil  table  prior  to  release.  Required  hammer 
heights  are  specified  in  MIL-S-901  based  on  the  total  weight 
attached  to  the  anvil  table. 

Deck  mounted  surface  ship  equipment  is  normally  attached 
to  the  anvil  table  using  the  standard  mounting  fixture  which 
is  designed  to  provide  a  great  deal  of  flexibility  in 
equipment  mounting  geometry.  The  standard  fixture  is  shown  in 
Figure  2.  Submarine  deck  mounted  equipment  is  attached  using 
the  coil  spring  soft  deck  simulator  shown  in  Figure  3.  When 
using  the  standard  fixture,  MIL-S-901D  specifies  the  number 
and  type  of  support  channels  to  be  used  based  on  the  equipment 
weight  and  mounting  bolt  spacing.  As  noted  by  Clements (1972) , 
the  arrangement  specified  in  MIL-S-901D  was  design  to  keep  the 
calculated  meiximum  stress  in  the  cheuinels  below  35,000  psi  in 
a  static  acceleration  field  of  50  g*s.  Although  not  by  design, 
using  this  arrangement  produces  an  equipment  excitation  in  the 
range  of  55  Hz  to  72  Hz.  The  soft  deck  simulator  has  a  natural 
frequency  in  the  reinge  of  20  Hz  to  25  Hz  and  is  designed  to 
simulate  the  natural  deck  frequencies  of  a  submarine. 

Equipment  tested  using  the  standard  fixture  will  be 
excited  predominately  at  a  single  frequency  in  the  55  Hz  to  72 
Hz  range  which  is  not  realistic  of  the  excitation  aboard  ship. 
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Figure  3.  Soft  Deck  Simulator  for  MWSM.  Courtesy  of  Hughes 
Aircraft . 
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Onboard  ship,  the  excitation  will  be  a  coit^lex  combination  of 
the  natural  frequencies  of  the  ship's  structure  in  the 
location  of  the  equipment  installation.  Shipboard  vibration 
may  excite  natural  frequencies  within  the  equipment  and  lead 
to  severe  damage  that  would  have  not  been  detected  on  the 
MWSM.  This  introduces  the  need  to  develop  a  replacement  for 
the  standard  fixture  which  will  produce  a  shock  excitation 
more  realistic  of  the  local  shipboard  response  to  underwater 
explosions . 

This  study  examines  the  design  of  a  multi-degree-of- 
freedom  fixture  to  be  installed  on  the  MWSM  which  can  be  tuned 
to  more  accurately  simulate  shipboard  shock  excitation  of 
combat  systems  equipment.  The  shock  spectrum  of  the  excitation 
input  to  three  representative  pieces  of  equipment  modeled  in 
the  DDG-51  Class  Ship  Pre-Shock  Trial  Analysis  performed  by 
Costanzo  and  Murray  (1991)  at  the  Underwater  Research  Division 
of  David  Taylor  Research  Center  (DTRC/UERD)  is  used  to 
identify  the  desired  excitation  response  of  the  fixture.  These 
shock  spectra  identify  the  frecjuency  components  of  the 
equipment  excitation  eind  provide  the  basis  for  tuning  the 
fixture.  This  fixture  is  proposed  as  an  economical  means  of 
providing  realistic  excitations  to  combat  systems  equipments 
when  shock  qpialifying  on  the  U.S.  Navy's  MWSM. 
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II.  BACKGROUND  PRESENTATION 


A,  THE  NEED  FOR  A  TUNED  MOUNTING  FIXTURE  FOR  THE  MWSM 

In  an  effort  to  the  increase  the  survivability  of  ships  in 
battle,  more  attention  is  being  focused  on  the  shock 
qualifications  of  combat  systems  equipment  installed  on  ship. 
To  ensure  the  reliability  of  this  equipment  in  a  mechanical 
shock  environment  it  must  first  be  "shock  qualified"  on  the 
appropriate  shock  machine.  For  the  test  to  be  sufficiently 
valid,  the  mechanical  shock  excitation  to  the  equipment  should 
be  simulated  as  closely  as  possible  to  the  expected  local  ship 
structure  response  to  an  underwater  explosion. 

MIL-S-901D  requires  that  the  test  item  be  mounted  to  the 
shock  machine  anvil  taible  in  a  manner  characteristic  of  its 
shipboard  orientation.  For  most  equipment,  this  requires  using 
a  specified  number  of  support  channels  and  rails  which  make  up 
the  standard  fixture  (Figure  2).  Chalmers  and  Shaw  (1989)  note 
that  most  users  of  MIL-S-901D  believe  that  using  the  standard 
fixture  on  the  M('7SM  produces  a  more  severe,  all  encompassing 
shock  excitation  than  would  actually  be  experienced  onboard 
ship.  They  provide  evidence  in  their  report  that  using  the 
standard  fixture  can  be  an  undertest  as  well  as  an  overtest 
since  a  high  frequency  mounting  will  indeed  pass  higher 
acceleration  levels  at  higher  frequencies,  but  will  not 
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provide  the  resonant  euiplif  ications  generated  by  a  lower 
frequency  fixture.  In  essence,  if  the  ship's  structure  excites 
the  equipment  at  a  lower  frequency  than  was  supplied  by  the 
MWSM,  the  installed  equipment  may  experience  severe  resonant 
vibrations  that  were  not  experienced  during  testing.  If  a  low 
frequency  fixture  such  as  the  soft  deck  simulator  were  used, 


it  too 

may  be  an  undertest 

since 

it  would 

not  excite 

the 

higher 

resonant  frequencies 

that 

may  be 

present  in 

the 

equipment . 

As  noted  in  Corbell ( 1992 ) ,  a  finite  element  transient 
shock  analysis  of  the  DDG-51  Class  Deck  House  was  conducted  by 
the  Ship  Structure  and  Protection  Department  of  David  Taylor 
Research  Center  (DTRC/UERD) .  The  preliminary  report  by 
Costanzo  and  Murray  (1991)  was  obtained  along  with  the 
predicted  shock  analysis  for  various  weight  combat  systems 
equipme**:  located  on  the  0-3  level  of  the  DDG-51  class  ship. 
From  this  report,  the  predicted  shock  excitations  to  three 
representative  combat  systems  equipments  were  used  as  the 
desired  response  for  a  tuned  fixture  on  the  MWSM.  The 
following  equipments  selected  fall  within  the  weight  range  of 
test  items  normally  qualified  on  the  MWSM: 

-  Radar  Receiver  Transmitter  (RT-1293 /SPS-67 )  325  lbs 

-  Beam  Programmer  {MX-10873/SPY-1D)  1000  lbs 

-  Radio  Frequency  Amplifier  (AM-7159/SPY-1B)  4600  lbs 
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The  predicted  acceleration  wave  forms  and  associated  shock 
spectra,  analyzed  for  maximum  shock  trial  severity,  for  the 
three  equipments  are  show  in  Figures  4  through  6 . 

The  shock  spectra  which  are  insensitive  to  small  waveform 
variations,  describe  the  characteristic  frequencies  of  the 
shock  induced  excitations.  The  shock  spectrum  is  generated  by 
plotting  the  maximum  absolute  response  of  a  single-DOF 
undamped  oscillator  as  a  function  of  its  natural  frequency 
when  subjected  to  the  base  excitation  of  interest.  It  is 
important  to  note,  that  when  attenpting  to  compare  rhock 
motion,  the  wave  "form"  in  the  time  domain  is  much  less 
important  than  the  wave  “characteristics"  displayed  in  the 
shock  spectrum  (frequency  domain).  Therefore,  when  trying  to 
reproduce  the  shipboard  excitation  on  the  MWSM  it  is  not 
necessary  to  reproduce  the  wave  form,  only  the  frequency 
components  at  proper  acceleration  amplitude  levels.  TVie 
predicted  shock  spectra  demonstrate  that  the  deck  house  shock 
environment  will  be  significantly  different  than  the  high 
energy,  single  frequency  shock  excitation  produced  by  using 
the  standard  fixture  on  the  MWSM.  MIL-S-901D  does  not  require 
that  any  particular  waveform  or  spectrum  be  reproduced, 
however  the  need  for  a  timed  mounting  fixture  for  the  MWSM 
clearly  exists . 
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Figure  5.  DTRC/UERD  Predicted  Beam  Programmer  Vertical 
Acceleration  Wave  Form  and  Shock  Spectrum. 
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PiSTure  6.  uT'-f'/'JERD  Predicted  Radio  Frequency  Amplifier 
Vertical  Acce  ’  e  .‘ation  Wave  Form  and  Shock  Spectrum. 
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B.  jmparing  shipboard  and  tuned  mounting  fixture  response 

Extending  the  concept  developed  by  Chalmers  and  Shaw 
(1989),  Corbell  proposed  that  a  two  degree-of-f reedom  (DOF) 
fixture  could  be  designed  for  the  MWSM  that  would  excite  an 
attached  piece  of  equipment  with  most  of  the  energy 
concentrated  at  two  resonant  frequencies.  If  the  fixture  was 
"tunable",  a  MIL-S-901D  user  could  refer  to  UNDEX  shoc)c 
spectra  data  or  modal  testing  data  for  a  specific  area  of  a 
ship  and  tune  the  fixture  to  simulate  two  of  the  dominant 
excitation  frequencies . 

Figure  7  shows  a  comparison  of  the  DTRC/UERD  DDG-51 
predicted  shocV:  spectrum  for  a  radar  receiver / transmitter  to 
the  computer  modeled  shoclc  spectra  of  the  Standard  Fixture  and 
two  DOF  fixture  on  the  MWSM.  The  DTRC/UERD  predictions  show 
significant  levels  of  acceleration  amplitudes  at  55  Hz  and  155 
Hz.  The  modeled  Standard  Fixture  provides  a  high  level  of 
acceleration  at  72  Hz  but  does  not  supply  a  sufficient 
amplitude  of  acceleration  between  130  Hz  and  170  Hz  to 
adequately  shoc)c  qualify  the  radar  with  respect  to  the 
DTRC/UERD  predictions. 

The  two  DOF  model  was  "tuned"  to  better  simulate  the  shoc)c 
characteristics  predicted  for  the  radar  receiver/ transmitter . 
As  can  be  seen  from  Figure  7,  the  two  DOF  model  reproduced  the 
equipment  accelerations  at  the  characteristic  frequencies  and 
cimplitudes  of  the  DTRC/UERD  predictions.  Equipment  qualified 
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Figure  7 .  Comparison  of  the  DTRC/UEPX)  Predicted  Vertical 
Acceleration  Shock  Spectra  to  the  Modeled  MWSM  Fixtures  for 
the  Radar  Receiver /Transmit ter  -  325  lbs 
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in  this  manner  could  be  highly  relied  on  to  perform  in  a 
shipboard  shock  environment. 

C.  TWO  DEGRE32--OF-FREEDOM  FIXTURE  MODEL 

1.  Coupled  and  Uncoupled  Natural  Frequencies 

The  Ml'JSM  with  equipment  attached  was  modeled  as  two 
mass -spring-damper  systems  coupled  together  as  shown  in  Figure 
8.  The  anvil,  when  struck  by  the  hammer,  experiences  a 
half-sine  wave  vertical  acceleration  impulse  of  approximately 
one  millisecond  in  duration.  The  magnitude  of  the  impulse  is 
controlled  by  the  hammer  height  prior  to  release.  This  impulse 
was  modeled  as  shown  in  Figure  9. 


TWO-DEGREE-OF-FREEDOM  FIXTURE  .MODEL 
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C2LiJ 
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K2 


Ml  =  Upper  Tier  Mass 
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Mounting  Hardware 
Effective  Spring  Mass 

M,  3=  Lower  Tier  Mass 
^  Includes:  Intermediate  Mass 

Effective  Spring  Mass 

K  =  Tier  Stiffness 
C  “  Tier  Strucrtural  Damping 
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Fisrure  8.  Modeled  Two  Degree  of  Freedom  Fixture  Subjected 
to  Base  Excitation. 
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MODELED  ANVIL  ACCELERATION  FOR  1  FOOT  HAMMER  DROP 


Figure  9.  Modeled  Anvil  Taible  Acceleration  for  a  One  Foot 
Hanuner  Drop . 


The 

mass  Ml 

represents 

the 

equipment  mass 

and 

associated 

mounting 

hardware 

while 

Mj  represent 

the 

intermediate  mass  2md  support  mountings .  Each  tier  possesses 
characteristic  stiffness  and  dairying  properties  designated  by 
K  cind  C  respectively. 

To  facilitate  the  solution  of  the  mathematical  model, 
the  equations  of  motion  for  the  system  can  be  expressed  in 
coordinates  relative  to  the  motion  of  the  base,  z.  Expressing 
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the  relative  coordinates,  and  yj  in  terms  of  absolute 

coordinates, 

yi=Xi-2 

(1) 

y,=Jc^-z 

y2=x2-z 

y-^^x^-t  (2) 

y2-^2-i 

yields  the  relative  coordinate  transformations  for 
displacement,  velocity  and  acceleration  respectively.  Solving 
for  the  coupled  equations  of  motion  gives  the  following  matrix 
equation : 


k  o] 

*  4* 

O 

y-\ 

-K 

(iCi+iCj) 


0 


By  asstiming  the  structural  damping  is  small,  the 
natural  frequencies  for  free  vibration  can  be  found  by 
assuming  the  damped  and  undanped  natural  frequencies  are 
approximately  equal.  Neglecting  the  deimping  matrix: 


k  o] 

k  -K 

H- 

o' 

o 

-kj^  {/Cj+iCj) 

0 

(4) 
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The  system  natural  frequencies  can  then  be  found  by  solving 
the  following  determinate  for  the  system  eigenvalues: 


det|[/C]  -CillM]  ]|=0 


where  o)^=2n/„ 


(5) 


Defining  the  ratio  between  the  upper  and  lower  tier 
masses  as. 


a  = 


5: 


(6) 


the  relationship  between  the  system's  coupled  and  uncoupled 
natural  frequencies^  can  be  described  by, 


V2 

4^=  —  v' A' f )  2 -  ( 2 f  1  2  (8) 


^  These  equations  are  based  on  the  valid  assumption 
that  the  structural  damping  is  small  and  that  the  damped  and 
undamped  natural  frequencies  are  equivalent . 
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where  f^i  and  fn2  are  the  coupled  system  natural  frequencies 
and  and  fj  are  the  uncoupled  tier  natural  frequencies 
defined  by: 


f  =~i_ 


5l 

M, 


^2  = 


2n\ 


h. 


(9) 


(10) 


If  the  desired  system  natural  frequencies  are  )cnown, 
equations  (7)  and  (8)  can  be  iteratively  solved  for  the 
required  uncoupled  natural  frequencies .  Knowing  the  mass 
associated  with  each  tier,  equations  (9)  and  (10)  can  be  used 
to  txine  the  fixture  to  the  required  stiffness,  K. 

From  the  shock  spectrum  shown  in  Figure  7,  Corbell 
chose  60  Hz  and  155  Hz  as  the  system  design  resonant 
frequencies  of  the  two  DOF  model.  Figure  10  represents  the 
iterative  solution  of  equations  (7)  and  (8)  for  the  uncoupled 
natural  frequencies  in  graphical  form.  As  an  example,  the  data 
was  plotted  as  a  function  of  the  mass  ratio  a=l  while  fixing 
the  lower  tier  natural  frequency  at  f2=100  Hz.  Choosing  an 
upper  tier  natural  frequency  of  94  Hz  from  the  coordinate  axis 
provides  the  desired  system  natural  frequencies  of  60  Hz  and 
155  Hz.  This  procedure  was  repeated  for  a  mass  ratio  of  0.7 
and  1.3  as  shown  in  Figures  11  and  12. 
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(ZH)  SHIDNanOHHH  TVmLVN  aHldilOD  WHISAS 


Figure  10 .  System  Coupled  Natural  Frequencies  as  a  Function 
of  the  Upper  Tier  Natural  Frequency  (fj  .  Lower  Tier  Natural 
Frequency  (fj)  Fixed  at  100  Hz.  ioass  Ratio  =  1. 
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UNCOUPLED  NATURAL  FREQUENCY  n  (HZ) 


^  in 
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(ZH)  SSIDKHnOHUd  TVmiVN  OHldAOD  KHISAS 

Figure  11.  System  Coupled  Natural  Frequencies  as  a  Function 
of  the  Upper  Tier  Natural  Frequency  (f^)  .  Lower  Tier  Natural 
Frequency  (f2)  Fixed  at  100  Hz.  Mass  Ratio  =  0.7. 
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(ZH)  SHONHQOHHd  TVHTLLVM  OHldTlOO  PMHJ.SAS 


Figure  12.  System  Coupled  Natural  Frequencies  as  a  Function 
of  the  Upper  Tier  Natural  Frequency  (f;)  .  Lower  Tier  Natural 
Frequency ( f 2 )  Fixed  at  100  Hz.  Mass  Ratio=1.3. 
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UNCOUPLED  NAIURAL  FREQUENCY  U  (HZ) 


2.  Analytical  Response  Calculations 

To  determine  the  response  of  the  two  DOF  fixture  to 
base  excitation,  the  matrix  equation  3  was  first  uncoupled 
using  the  normalized  transformation  matrix, 


(11) 


where  the  columns  of  [O]  are  the  mode  shape  vectors 
(eigenvectors)  .  The  resulting  uncoupled  equations  of  motion  in 
natural  coordinates  are: 

fii+2Coj^jf|i+o)LTi^  =  -Zi  (12) 

ti2+2Co)„2f|2+wLT)2  =  --22  <13  ) 

where  ^  is  the  critical  damping  ratio  estimated  from  the  MWSM 
calibration  test  data  cortpiled  by  Costanzo  and  Clements 
(1988) .  Examination  of  the  calibration  data  showed  that  the 
test  weight  acceleration  response  damped  out  in  approximately 
.3  to  .5  seconds  after  hammer  iirpact.  This  corresponds  to  a 
critical  dairying  ratio  of  between  three  to  five  percent. 
Figures  13  through  15  show  the  damped  acceleration  response 
for  the  two  DOF  model  for  three  values  of  the  critical  danping 
ratio . 
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Figure  13 .  Modeled  Two  DOF  Fixture  Acceleration  Response 
for  Zeta  =  0.03. 
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TIME  (sec) 


Flgxirs  15.  Modeled  Two  DOF  Fixture  Acceleration  Response  for 
Zeta  =5  0.05. 


TIME  (sec) 


Equations  (12)  and  (13)  where  solved  independently 
using  a  continuous-time  linear  system  simulation.  The  solution 
in  natural  coordinates  was  then  transformed  into  relative 
coordinates  using  equation  (11)  and  then  to  absolute 
coordinates  using  equation  (1)  and  (2)  to  obtain  the  system 
displacement,  velocity  and  acceleration  responses.  The  MATLAB™ 
algorithm  solving  the  two  DOF  system  is  presented  in 
Appendix  A. 

The  two  DOF  model  and  numerical  response  algorithm 
provided  the  prerequisite  information  needed  for  designing  and 
predicting  the  response  of  the  proposed  fixture.  In  Chapter 
IV,  a  semi-definite  three  degree  of  freedom  model  will  be 
introduced  which  accounts  for  the  unconstrained  4500  lb  anvil 
table  interaction  with  the  fixture  after  the  time  of  hammer 
impact . 
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III.  TWO  DEGREE -OF “FREEDOM  TUNED  FIXTURE  DESIGN 

The  designs  discussed  in  this  chapter  are  based  on  the  two 
DOF  model  presented  in  Chapter  II.  Tne  procedure  followed  vas 
to  first  identify  the  desired  system  natural  frequencies  for 
the  fixture  and  then  determine  the  required  tier  natur>^x 
frequencies  using  equations  (7)  and  (8)  .  Once  the  tier  natural 
frequencies  were  known,  the  required  stiffness  for  each  tier 
could  be  determined  using  equations  (9)  and  (10).  The  mass, 
stiffness  and  damping  for  each  tier  v;ere  input  into  a 
niimerical  algorithm  to  determine  the  fixture  response  to  the 
half-sine  wave  impulse  shown  in  Figure  9,  The  displacements  of 
each  tier  relative  to  the  anvil  table  were  determined  for  a 
variety  of  hammer  heights,  equipment  weights  and  system 
natural  frequencies.  The  maximum  dynamic  relative  displacement 
within  the  fixture  was  obtained  for  the  above  conditions  and 
then  used  to  conduct  a  static  stress  analysis  of  the  fixture. 
A  detailed  procedure  of  the  design  process  is  presented  below. 

A.  DESIGN  OF  THE  TWO  DOF  TUNED  FIXTURE 


1.  Design  Considerations 


The 

des 

ign 

phase  of 

the 

research  began 

with 

identifying 

the 

two 

resonant 

frequencies  desired  for 

the 

fixture.  After 

referring  with 

Mark 

McClean  at  Naval 

Sea 

Systems  Command,  the  previously  desired  frequencies  of  60  Hz 
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and  155  Hz  had  been  lowered  to  30  Hz  and  80  Hz  for  the  goal 
frequencies.  Modal  testing  in  the  applicable  areas  of  the  ship 
indicated  that  the  shipboard  frequencies  were  lower  than  the 
DTRC/UERD  predictions. 

A  second  design  consideration  was  cost.  The  Navy  could 
easily  develop  a  pulse-shape  matching  machine  capable  of 
reproducing  shipboard  shock  spectra  but  the  cost  per  machine 
would  be  significant.  The  Navy  currently  has  13  MWSM  in  use 
and  replacing  them  with  million  dollar  machines  is  not 
economically  feasible.  The  design  approach  was  to  develop  a 
fixture  that  used  many  of  the  current  parts  and  accessories 
already  in  use  on  the  MWSM.  By  doing  so,  fabrication  and 
material  cost  could  be  minimized. 

MIL-S-901D  provides  some  general  guidelines  for 
designers  building  equipment  to  endure  shock  excitations,  many 
of  which  are  applicable  to  the  fixture  design.  Materials 
recommended  for  structural  members  are  those  with  high  yield 
strength,  high  ductility,  high  fracture  touglmess  and  when 
possible,  light  weight.  Cast  iron  and  cast  aluminum  have 
generally  proven  to  be  unsatisfactory  when  used  as  strength 
members  due  to  their  high  notch  sensitivity  cind  brittleness. 
Areas  normally  subjected  to  conpressive  stress  under  static 
loading  conditions  may  experience  tensile  stress  when 
subjective  to  the  cyclic  loading  experienced  under  shock  and 
vibration.  This  fact  must  be  kept  in  mind  specifically  when 
designing  welded  and  bolted  joints.  These  joints  tend  to  cause 
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stress  concentrations  which  may  be  acceptable  in  compression 
but  may  lead  to  failure  in  tension. 

2.  Preliminary  Design 

The  following  designs  were  based  on  the  two  degree  of 
freedom  model  shown  in  Figure  8.  The  focus  here  was  to  define 
a  structure  with  a  suitable  set  of  tier  stiffness  and  mass 
components  which  would  provide  the  desired  response  without 
yielding  or  structural  failure.  Numerous  designs  were  explored 
and  then  narrowed  to  two  for  this  research. 
a.  Coil  Spring  Deck  Deelga 

The  soft  deck  simulator  shown  in  Figure  3  provides 
some  obvious  quali.ties  which  would  be  quite  desirable  for  this 
design.  The  most  significant  quality  being  the  ease  in  tuning 
the  fixture.  As  the  weight  of  the  attached  equipment  varies 
the  fixture  is  easily  tuned  by  adding  or  removing  spring 
cartridges  to  maintain  the  desired  natural  frequency.  A  second 
quality  is  the  relatively  large  vertical  displacement  allowed 
within  the  fixture.  Since  the  relative  displacement  is 
inversely  proportional  to  square  of  the  natural  frequency,  it 
becomes  increasingly  necessary  to  allow  more  relative  motion 
between  the  tiers  as  lower  frequencies  are  trying  to  be 
obtained.  The  coil  springs  will  allow  up  to  two  inches  of 
travel  without  being  over  stressed  or  becoming  coil  bound. 

Another  feature  of  using  coil  springs  is  the 
ability  to  distribute  the  spring  force  over  the  entire  surface 
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area  of  the  table  top  which  has  two  advantages.  First,  as 
mentioned  in  Clements  (15)72),  the  MWSM  is  not  a  perfect 
machine  and  does  not  alvays  produce  a  strictly  vertical 
impulse  on  the  anvil  table.  By  distributing  the  spring  force, 
the  fixture  becomes  more  stable  since  it  can  better  resist 
the  rotational  moment  produced  by  the  supported  mass  when  its 
center  of  gravity  is  not  in  line  with  the  direction  of  the 
impulse  force.  Secondly,  since  the  force  is  distributed,  the 
bending  moments  applied  to  the  rigid  intermediate  mass  (Mj) 
are  minimized.  This  reduces  unwanted  secondary  frequencies 
generated  by  deflections  within  the  structure  of  the 
intermediate  mass. 

To  achieve  the  desired  system  natural  frequencies 
of  30  Hz  and  30  Hz  the  required  uncoupled  frequencies  from 
equations  (7)  and  (8)  were  determined  to  be  fi=48  Hz  and  f2=50 
Hz .  The  mass  ratio  a  was  set  equal  to  one  with  each  tier  mass 
weighing  3000  lbs.  This  would  leave  1400  lbs  available  of  the 
7400  lbs  allowable  on  the  MWSM  for  mounting  hardware  and  the 
inclining  fixture  called  for  in  MIL-S-901D.  Using  equations 
(9)  and  (10)  the  required  upper  and  lower  tier  stiffness  was 
determined  to  be  Kj=7 . 66x10®  Ib/in  and  K2=7. 06x10®  Ib/in. 

The  coil  springs  currently  used  in  the  soft  deck 
simulator  have  a  stiffness*  of  2205  Ib/in  and  weigh 

*  This  is  the  mean  stiffness  per  coil  spring  as 
determined  by  Steve  Schecter  at  Hughes  Air  Craft,  Fullerton, 
Ca . 


approximately  15.5  lbs  each.  This  would  require  348  springs  on 
the  upper  tier  and  320  springs  on  the  lower  tier  for  a 

combined  weight  of  10,354  lbs  for  the  springs  alone.  This  is 

obviously  an  unacceptable  weight  for  use  on  the  MWSM. 

An  alternate  off-the-shelf  spring  was  selected  with 
a  stiffness  of  6090  Ib/in,  a  weight  of  18  lb  per  spring  and 
similar  dimensions.  This  would  require  126  springs  on  the 
upper  tier  and  116  springs  on  the  lower  tier  for  a  total 
weight  of  3280  lbs.  This  reduced  the  spring  weight 
significantly  but  still  consvuned  too  much  of  the  weight 

capacity  of  the  MWSM.  The  advantages  of  the  coil  spring  design 
had  to  be  abandoned  for  a  design  with  a  higher 

stiff ness-to-weight  spring  mechanism.  The  beeun  loading  concept 
used  in  the  standard  fixture  was  a  logical  choice. 
b.  Beam  Spring  Design 

■  Obtaining  the  necessary  stiffness  while  minimizing 
weight  encouraged  the  use  of  beams  for  the  spring  elements. 
The  beam  elements  were  modeled  as  being  siimly  supported  with 
symmetric  loading  as  shown  in  Figure  16.  The  stiffness  of  the 
beam  is  given  by  equation  (14) . 


X=£  =  _lil_ 

y  a2{3L-4a) 


(14) 
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Simple  Supports  -  Symetrlc  Loads 


Figure  16 .  Beam  Bending  Model  for  Determining  Beam 
Stiffness  . 

This  design  was  based  on  using  the  existing 
shipbuilding  rails  as  the  spring  beam  supports .  From  the 
drawings  of  the  MWSM  in  MIL--S-901D,  this  fixed  L  at 
approximately  50  inches.  The  load  spacing  was  set  to  24  inches 
making  a=13  inches.  Young's  Modules  was  taken  to  be  E*30xl0‘ 
psi  for  steel  and  the  area  moment  of  inertia  was  set  to  9.2 
in*  ,  based  on  the  current  back-to-back  C  channels  (4x7.25#) 
being  used  on  the  standard  fixture.  The  result  was  a  stiffness 
of  2 .17x10®  lb/ in  with  each  beam  set  weighing  approximately  70 
lb.  To  achieve  the  desired  natural  frequencies,  the  lower  tier 
required  4  beam  sets  ^u:ld  the  upper  tier  required  3  beam  sets 
for  a  total  weight  of  490  lbs.  The  beams  provided  a 
significeuit  weight  savings  and  became  the  focus  of  the 
detailed  design. 


B.  DETAIL  DESIGN  OF  THE  BEAM  SPRING  FIXTURT^ 

In  designing  the  two  DOF  fixture,  priority  was  placed  on 
the  ability  to  tune  the  fixture  to  achieve  the  desired  system 
natural  frequencies  over  the  widest  possible  range  of 
equipment  weights.  An  equally  important  consideration  was 
minimizing  the  weight  of  the  fixture  while  still  maintaining 
the  characteristics  of  the  shipboard  shock  wave.  A  1/4  scale 
model  of  the  proposed  two  DOF  fixture  is  presented  in  Figure 
17.  The  f>  llowing  paragraphs  highlight  the  significant  factors 
dealt  with  during  the  design  process. 

1 .  Choice  of  Beams 

The  back-to-back  C  chainnels  used  or  the  standard 
fixture  were  first  investigated  for  use  in  the  two  DOF  design. 
These  channels  would  mount  to  the  shipbuilding  support  rails 
in  the  current  fashion  using  the  existing  claitping  method. 
However,  in  an  effort  to  simplify  the  design  and  to  use  a 
symmetric  cross  section,  I-beams  were  ultimately  chosen.  The 
I-beams  used  in  the  design  are  W4xl3  beams  with  an  area  moment 
of  inertia  of  11.3  in*  resulting  in  a  stiffness  of  2.46x10* 
lb/ in.  By  modifying  the  existing  clcunping  design,  the  fixture 
could  be  assembled  without  requiring  emy  drilling  or  welding 
of  the  I  beams  in  areas  of  high  stress.  To  achieve  the  system 
natural  frequencies  desired,  required  four  I-beams  on  the 
bottom  tier  and  three  on  the  top  tier.  To  add  more  flexibility 
in  mounting  the  equipment  to  the  fixture,  the  design  allows 
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for  interchanging  the  I-beams  with  the  back-to-back  C  channels 
on  the  upper  tier . 

2 .  Tuning  The  Fixture 

As  the  equipment  weight  attached  to  the  fixture 
changes,  the  fixture  is  tuned  by  altering  the 
stif fness-to-mass  ratio  of  the  upper  euid  lower  tiers.  On  the 
upper  tier,  the  stiffness  is  increased  by  reducing  the 
dimension  "a“  shown  in  Figure  16.  The  initial  design  of  the 
lower  tier  allowed  for  the  same  adjustment  to  alter  the 
stiffness.  In  an  effort  to  increase  the  rigidity  of  the 
intermediate  mass  Mj,  and  to  increase  stability  by  lowering 
the  center  of  gravity  of  the  entire  fixture,  the  design  was 
changed  fixing  a=13.5  inches  for  the  lower  tier  beams.  The 
lower  tier  natural  frequency  is  adjusted  by  adding  or  removing 
ballast  to  the  intermediate  mass  as  necessary. 

When  tuning  the  fixture,  the  effective  mass  of  the 
springs  was  taken  into  account.  Using  Furtis  (1972)  as  a 
reference,  the  effective  mass  of  the  spring  elements  was 
determined  to  be  52%  of  the  total  beam  mass.  For  Kj,  this  mass 
was  added  to  the  intermediate  mass  Mj.  For  Kj,  52%  of  the 
spring  mass  was  added  to  the  upper  tier  mass  and  the 
remaining  48%  added  to  Mj  since  it  would  be  in  motion  with  the 
intermediate  mass.  The  distribution  of  the  effective  spring 
mass  is  some  what  trivial  since  its  magnitude  is  quite  small 
relative  to  the  tier  masses. 


3.  Intermediate  Mass  (M,)  Design 

As  shown  earlier  in  equations  (7)  and  (8),  the  ratio 
of  the  upper  tier  mass  to  lower  tier  mass  (a)  has  an  effect  on 
the  system  response  to  the  uncoupled  tier  frequencies.  The 
plot  in  Figure  18  depicts  this  relationship.  For  desired 
system  natural  frequencies  of  30  Hz  and  80  Hz,  the  plot  shows 
the  required  tier  natural  frequencies  as  a  function  of  a.  As 
a  increases,  the  required  upper  tier  natural  frequency 
decreases  while  the  lower  tier  increases .  For  a  greater  than 
1.3,  the  required  uncoupled  natural  frequencies  diverge  very 
quickly  indicating  that  the  desired  system  response  is  not 
achievable.  This  signifies  that  the  intermediate  mass,  Mj, 
must  be  greater  than  or  equal  to  77%  of  the  equipment  mass. 
The  design  weight  of  Mj  including  the  mounting  hardware  is 
17  00  lbs  and  will  require  adding  ballast  for  equipment  weights 
over  2200  lbs.  The  result  is,  Mj  will  consume  a  significant 
portion  of  the  weight  capacity  of  the  MWSM  and  will  limit 
equipment  weight  to  approximately  3500  lbs.  This  weight 
allowance,  although  somewhat  reduced,  will  still  allow 
significant  numbers  of  equipments  to  be  qualified  with  this 
design . 

A  second  design  consideration  for  M,  was  that  it  must 
be  stiff  enough  to  appear  as  a  rigid  body  to  the  rest  of  the 
fixture.  This  would  prevent  unwanted  deflections  and 
vibrations  not  associated  with  the  desired  system  response. 
The  necessary  stiffness  was  achieved  by  forming  a  grid 
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Figure  18.  Required  Upper  and  Lower  Tier  Natioxal  Frequencies 
as  a  Funct:ion  o£  the  Hass  Ratio. 
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Mass  Ratio  (Mi/M2) 


structure  of  I-beams  with  an  overall  stiffness  of  l.OxlO’ 
Ib/in  which  is  approximately  13  times  higher  than  Kj  or  Kj . 

To  preclude  the  need  for  fabricating  a  different 
spring  beam  for  the  upper  tier  and  lower  tier,  the 
intermediate  mass  provides  the  same  mounting  geometry  as 
the  shipbuilding  support  rails  mounted  on  the  MWSM  anvil 
table.  This  allows  interchanging  the  upper  and  lower  tier 
springs  as  well  as  the  back-to-back  C  channels  currently  used 
as  part  of  the  standard  fixture.  In  some  instances,  mounting 
equipment  on  the  upper  tier  with  back-to-back  C  channels  will 
be  more  convenient  and  provide  more  flexibility  in  cuning  the 
fixture . 

C.  STRESS  ANALYSIS  OF  THE  TWO  DOF  FIXTURE 

Considering  only  deflections  in  the  vertical  direction, 
evaluation  of  the  most  likely  modes  of  failure  focused  on  the 
tensile  stress  in  the  fasteners  and  bending  stress  in  the 
spring  beams  and  clamps.  The  numeric  model  was  loaded  to  its 
weight  capacity  and  then  excited  with  the  maximum  expected 
acceleration  level.  The  upper  tier  weight  was  set  to  3500  lbs 
and  the  lower  tier  weight  set  to  2700  lbs.  Using  the  300  Hz 
filtered  data  shown  in  Figure  19,  the  maximum  expected  anvil 
acceleration  was  picked  to  be  200  g's  for  a  five  foot  hammer 
drop.  Using  these  parameters  in  the  numeric  model,  the 
relative  displacements  of  the  upper  and  lower  tiers  with 
respect  to  the  anvil  table  were  determined.  Figure  20  is  a 
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(6)  N0ilVb3'!300'7  3ieVi-HANV  XV3d 


Figure  19.  MWSM  Peak  Anvil  Table  Acceleration  vs.  Hammer 
Height.  Courtesy  of  Clemenets  (1972). 
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HAMMER  IMPACT  VELOCIT 


spring  beam  can  be  calculated. 


Af=Fa 


(16) 


Knowing  the  bending  moment  ,M,  the  maximum  bending  stress 
was  determined  from  equation  (17); 


o 


t>ending' 


Me 

I 


(17) 


where  c=2 . 0  inches  for  the  four  inch  I-beam.  The  shear  stress 
in  tne  beams  was  determined  to  be  negligible  and  not  included 
in  the  calculations.  (Shear  stress  is  zero  at  the  extreme 
fiber  were  the  bending  stress  is  maximum) . 
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The  equivalent  static  force  ,F,  was  used  correspondingly 
to  determine  the  stresses  in  the  clanps  and  bolts  fastening 
the  fixture  together.  The  maximiom  stress  in  the  spring  beams 
was  determined  to  be  47  kpsi  and  13  kpsi  in  the  clamps. 
Maximum  tensile  stress  in  the  bolts  was  determined  to  be  7.5 
kpsi  . 

D.  DESIGN  DRAWINGS  AND  MATERIAL  SELECTION 

Detailed  design  CAD  drawings  of  the  proposed  two  DOF  tuned 
fixture  were  provided  to  Naval  Sea  Systems  Command  for 
approval  and  fabrication.  A  set  of  these  drawings  have  been 
included  in  Appendix  B  of  this  document. 

The  material  selected  for  the  spring  beams  has  a  yield 
strength  of  50  kpsi.  Materials  selected  for  all  other 
structural  elements  is  standard  36  kpsi  yielc*  strength  mild 
steel.  Bolt  material  is  high  quality  grade  5  carbon  steel 
alloy.  All  aaterials  have  been  received  and  fabrication 
expected  to  begin  by  1  June  1993.  Upon  completion,  the  fixture 
will  be  sent  to  Naval  Underwater  Systems  Center,  New  London, 
Ct .  for  testing. 

E.  PREDICTED  RESPONSE  FOR  THE  TWO  DOF  FIXnjRE 

To  predict  the  performance  of  the  proposed  two  DOF 
fixture,  a  sensitivity  analysis  of  the  model  was  performed  to 
determine  the  tunability  and  the  response  characteristics  of 
the  fixture  as  equipment  weight  was  varied.  For  the  analysis. 
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the  desired  system  natural  frequencies  were  chosen  to  be  20  to 
40  Hz  for  fnj  2ind  7  0  to  30  Hz  for  f^.  The  equipment  weight 
including  the  mounting  hardware  varied  from  500  lb  to  3500  lb. 

1.  Tunability  of  the  Fixture 

To  examine  the  txinability  of  the  fixture,  a  worksheet 
was  produced  in  Mathcad®  which  iteratively  solves  equations 
(7)  and  (8)  for  the  required  uncoupled  tier  frequencies  f^  and 
fj  given  the  desired  system  natural  frequencies  f„i  amd  f,^2  and 
the  mass  ratio  a.  Using  the  worksheet,  the  required  tuning  can 
easily  be  determined  as  equipment  weight  or  desired  system 
natural  frequencies  change.  A  copy  of  this  worksheet  is 
provided  in  Appendix  C.  Figure  21  shows  the  achievcible  system 
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Figure  21.  System  Natural  Frequencies  as  a  Function  of  the 
Upper  Tier  Weight 


TABLE  I .  TUNING  FOR  VARIOUS  UPPER  TIER  WEIGHTS  . 


2-Upper  tier  equiment  mounting  spacing.  Refer  to  ficpare  16. 


natural  frequencies  as  the  weight  of  the  upper  tier  is  varied 
between  500  and  3  500  lbs.  TABLE  I  shows  the  necessary 
adjustments  which  would  be  required  to  tune  the  fixture.  The 
sharp  rise  in  system  natural  frequencies  for  upper  tier 
weights  less  thaui  1500  lbs  indicates  that  adding  ballast  to 
the  upper  tier  may  be  required  for  some  light  weight 
equipment . 
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A  second  analysis  was  performed  to  determine  the  range  of 
system  natural  frequencies  obtaineOole  for  a  specified 
equipment  weight.  The  equipment  weight  was  given  to  be  2500 
lbs  and  the  fixture  tuned  by  varying  the  weight  of  Mj  and 
adjusting  the  equipment  mounting  spacing  ("a"  in  Figure  16). 
For  an  equipment  weight  of  2500  lbs.  Figure  22  demonstrates 
the  obtaincible  system  natural  frequencies.  The  data  provided 
in  TABLE  II  indicates  the  necessary  adjustments.  Figures  21 
and  22  demonstrate  that  this  fixture  design  could  be  tuned  to 
the  desired  system  natural  frequencies  under  a  variety  of 
loading  conditions. 
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TABLE  II.  TUNING  FOR  AN  UPPER  TIER  WEIGHT  OF  2500  LBS. 


Case 

Number 

Lower  Tier 
Weight 
(lb) 

•’a"  Upper 
Tier 
( inches) 

#  Beams 
Upper 
Tier 

fn, 

(Hz) 

fnj 

(HZ) 

1 

3000 

16 

2 

31 

70 

2 

3000 

15 

2 

31 

71 

3 

2750 

16 

2 

31 

73 

4 

2750 

15 

2 

32 

75 

5 

2500 

16 

2 

31 

76 

6 

2250 

15 

2 

32 

81 

7 

2250 

14 

2 

33 

83 

8 

2000 

13 

2 

34 

89 

9 

2000 

14 

3 

37 

96 

10 

1750 

3 

38 

107 
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2.  Two  DOF  Fixture  Response  Characteristics 

As  mentioned  in  Chapter  II,  for  the  fixture  to  be 
effective,  it  should  provide  the  equipment  with  accelerations 
that  are  "characteristic"  of  the  expected  shipboard 
excitations.  Since  the  desired  frequencies  have  been  lowered 
by  NSWC  to  3  0  Hz  cind  80  Hz,  it  is  difficult  to  compare  the 
modeled  results  to  the  DTRC/UERD  predictions.  However,  as 
shown  for  five  different  equipment  weights.  Figures  23  through 
2  8  demonstrate  that  the  modeled  time  history  accelerations  and 
associated  shock  spectra  provide  significant  amplitudes  of 
acceleration  at  the  desired  frequencies. 
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ACCEI.EKATION  (g’s)  ACCELERATION  (g's) 


MODELED  TWO  DOF  FDCTURE  UPPER  TIER  ACCELERATION 


modeled  TWO  DOF  FIXTURE  ACCELERATION  SHOCK  SPECTRA 

450, - - - - - 


Figure  23.  Predicted  Acceleration  Wave  Fonti  and  Shock 


Spectrum  for  1500  lb  Upper  Tier  Weight. 
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ACCELERATION  (g‘s)  ACCELERATION  (g’s) 


MODELED  TWO  DOF  FDCTURE  UPPER  TIER  ACCELERATION 


MODELED  TWO  DOF  FIXTURE  ACCELERATION  SHOCK  SPECTRA 


FREQUENCY  (Hz) 


Figure  24.  Predicted  Acceleration  Wave  Form  and  Shock 
Spectmm  for  2000  lb  Upper  Tier  Weight. 
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ACCELERATION  (g's)  ACCEI.ERATION  (g's) 


Spectrtim  for  2500  lb  Upper  Tier  Weight. 
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ACCELERATION  (g's)  ACCELERATION  (g’s) 


Spectrum  for  3000  lb  Upper  Tier  Weight. 
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ACX:ELERATJ0N  (g’s)  ACCEI.ERATJON  (g's) 


MODELED  TWO  DOF  FIXTURE  UPPER  TIER  ACCELERATION 


5  Foot  Hammer  Drop  (200  g's) 
_  Upper  Tier  Weight  =  3500  Ib 
Lo^r  Tier  Weight  =  2700  lb 
Zeta  =  0.04 
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Figure  27 .  Predicted  Acceleration  Wave  Form  and  Shock 
Spectrum  for  3500  lb  Upper  Tier  Weight. 
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IV.  SEMI -DEFINITE  THREE  DEGREE -OF -FREEDOM  MODEL 

As  mentioned  in  Chapter  II,  Corbell  formulated  a  two- 
degree-of- freedom  model  to  simulate  the  fixture  response  to 
the  anvil  table  input.  The  fixture  model  was  shown  in  Figure 
8.  The  input  into  the  model  was  the  half-sine  wave  base 
acceleration  impulse  described  in  Figure  9.  After  the  base 
impulse  passed  {1  msec),  the  anvil  table  was  considered 
stationary  and  the  fixture  vibrated  as  a  two  DOF  mass-spring- 
demiper  system  attached  to  a  fixed  foundation.  This  model 
provided  useful  information  required  for  determining  the 
feasibility  and  design  of  the  proposed  fixture.  In  this 
chapter,  a  semi-definite  three  DOF  model  will  be  introduced 
which  will  take  into  account,  the  4500  lb  anvil  table 
interaction  with  the  fixture  after  the  time  of  hammer  impulse. 

A.  SEMI -DEFINITE  THREE  DOF  MODEL  DERIVATION 
1.  Model  Description 

As  presented  in  Clements  (1972),  the  anvil  table  is 
bolted  to  the  machine  foundation  with  12,  2-inch-diameter 
bolts  in  a  manner  which  allows  the  table  three  inches  of 
vertical  travel  after  hammer  impact.  This  vertical  distance 
can  be  decreased  by  raising  the  table  with  pneumatic  jacks 
prior  to  performing  the  test.  Approximately  50  msec  after 
hammer  impact,  the  table  reaches  the  limit  of  vertical  travel 
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and  experiences  a  sharp  "table  reversal"  in  the  form  of  a 
fairly  simple  negative  acceleration  pulse.  Since  the  duration 
of  the  reversal  is  2  to  4  times  that  of  the  hammer  impact  and 
energy  has  been  expended  ir  the  system,  the  magnitude  of  the 
acceleration  pulse  is  much  less  than  that  caused  by  the 
hammer.  As  the  table  falls  back  down  on  the  foundation  it 
experiences  another  impulse  of  longer  duration  and 
consequently  lower  acceleration.  From  the  time  of  hammer 
impact  to  the  time  of  table  reversal,  the  anvil  table  is 
unconstrained  and  the  entire  structure  behaves  as  the  semi- 
definite  three  DOF  system  shown  in  Figure  28.  Considering  the 
initial  impact  as  significantly  more  sever  than  the 
subsequent,  it  was  possible  to  simplify  the  model  and  not 
include  the  table  stop  in  the  mathematical  analysis  for 
determining  the  characteristics  of  the  fixture  response. 

Referring  to  Figure  28,  the  "upper  tier"  mass  Mj 
represents  the  mass  of  the  equipment  and  mounting  hardware  and 
the  effective  spring  mass  for  the  upper  tier.  represents 
the  intermediate  mass  and  the  effective  spring  mass  for  the 
lower  tier.  Mj  is  the  mass  of  the  anvil  table.  The  spring 
stiffness  is  designated  by  Kj  for  the  upper  tier  and  Kj  for  the 
lower.  Again,  as  in  the  two  DOF  model,  the  structural  damping 
is  considered  proportional  to  the  mass  and  stiffness.  The 
upper  tier  and  lower  tier  damping  is  designated  by  Cj  and  Cj 
respective  .'.y . 
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THREE  DOF  MODEL 


Three  DOF  Model  for  Fixture  on  MWSM 


Referring  to  Clements  (1972),  the  peak  anvil  table 
acceleration  produced  by  the  half-sine  impulse  is  very  nearly 
a  linear  function  of  the  hammer  impact  velocity  and  is 
essentially  "independent"  of  the  load  when  it  is  channel 
mounted.  Using  this  information,  the  forcing  function  F(t)  can 
be  modeled  as  simply  the  anvil  acceleration  times  the  anvil 
mass  and  is  independent  of  the  load  attached. 


Fit)  (18) 

Since  the  test  scnedule  for  the  MW3M  listed  in  MIL-S- 
901D  specifies  the  "hammer  drop  height"  for  a  given  table 
weight,  it  is  convenient  to  write  F(t)  as  a  function  of  the 
height  of  the  hammer  drop.  Referring  back  to  Figure  19,  the 
300  CPS  filtered  data  shows  a  linear  correlation  between  anvil 
acceleration  and  velocity  with  a  slope  of  approximately  11.35 
g' sec/ft.  Using  the  relation,  V=V2gh,  the  peak  anvil  table 
acceleration  can  be  related  to  the  hammer  height  by  the 
followii  relation: 


Ji^ipeak)  =(11.35)  (32.2)/?^ 
where:  h=Haininer  Height  Above  Anvil  Table  in  Feet 


(19) 
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Knowing  the  peak  acceleration  of  the  half  sine  wave 


impulse,  the  forcing  function  F(t)  can  then  be  modeled  as; 


F(C]  =293  3v^(Ar3)  sin(G)t)  0;Jti0.001  sec 

F(t)-0  t>0.001  sec  (20) 

where:  <.)=l000*r 

The  1000  CPS  filtered  data  shown  in  Figure  19 
represents  the  table’s  750  Hz  longitudinal  mode  of  damped 
vibration  which  last  for  only  about  5  periods  and  does  not 
significantly  contribute  to  the  motion  of  the  fixture. 

2.  Mathematical  Analysis  of  the  Three  DOF  Model 

The  number  of  degrees  of  freedom  of  a  system  is 
defined  as  the  number  of  independent  coordinates  necessary  to 
describe  the  motion  of  a  system  completely.  For  the  model 
shown  in  Figure  28,  this  will  recruire  three  independent 
coordinates  to  define  the  motion  of  the  three  masses  of  the 
system.  These  coordinates,  Xj,  and  Xj,  define  the  absolute 
displacements  of  the  masses,  Mj,  Mj  and  Mj,  respectively.  The 
time  derivatives  of  these  coordinates  yield  the  velocity  and 
acceleration: 


1^2  '  ^2 

(Displacement) 

,X2  ,X^ 

( Veloci  ty) 

(21) 

1  X2  1  ^3 

(Accel  elation) 
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The  resulting  equations  of  motion  in  absolute 
coordinates  can  be  written  as : 


^kj^x^ -c^X2  =0  ( 22 ) 

1112X2  +  {C2+C2)  X2+  {k^+k2)  X2- c^x^ - C2X3 -/C2X3 =0  (23) 

^3X3+03X3  +^2-^3  "^2^2  "■^2-^2  (  24  } 

The  above  equations  of  motion  can  be  conveniently 
expressed  in  the  following  matrix  form: 


■^1 

0 

O' 

O' 

f'] 

O' 

0] 

0 

0 

^2 

.  + 

-^1 

(C3+C2) 

-<^2 

.  -f 

-K 

{k,^k2) 

-^2 

^2 

0 

0 

0 

m3 

^3, 

0 

-^2 

•^2 

0 

-^2 

K. 

1:^3 

Since  the  structural  damping  of  the  fixture  is  small, 
the  system  natural  frequencies  (free  vibration)  can  be  found 
by  assuming  that  the  damping  is  negligible  and  that  the  damped 
and  undamped  frequencies  are  equal.  Neglecting  danping  and 
assximing  free  vibration,  equation  (25)  can  be  rewritten  as: 


0  O' 

■  k^  -k^  0 

o' 

0 

0 

*2 

-k^  ( /Cj  +ic2 )  -k2 

X2 

= 

0 

^  0  0  m3 

.^3, 

0  -iCj  k2 

;^3 

0. 

(26) 
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The  system  natural  frequencies  are  then  found  by 
solving  the  following  determinant  for  the  system  eigenvalues: 


det|[/0  ]|  =  0 

where  (j3„=2nf„ 


(27) 


The  solution  to  equation  (27)  is  shown  below  in 
equations  (28)  through  (30)  revealing  that  one  of  the  three 
system  natural  frequencies  is  zero.  This  zero  frequency 
corresponds  to  che  rigid  body  mode  of  vibration. 


{Rigid  Body  Mode) 


(28) 


0} 


n2' 


1 

B-* 

^|  2 

N 

(29) 


1 

B-*-. 

2 

N 

(30) 


where 


M,"  M, 
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As  described  in  Chapter  II,  if  the  desired  system 
coupled  natural  frequencies  fnj.f^j  and  f^j  are  known,  equations 
(29)  and  (30)  can  be  iteratively  solved  to  determine  the 
required  uncoupled  tier  natural  frequencies  f i  and  f^ .  A 
Mathcad®  program  is  listed  in  Appendix  C  which,  given  the 
c’-'sired  system  natural  frequencies  v/ill  solve  for  the  required 
tier  frequencies  . 

Having  found  the  system  natural  frequencies,  the 
natural  modes  of  vibration  (eigenvectors)  can  be  solved  for 
using  equation  (31). 

-^1  C 

-iCj  (ki+kj) -k2 
0  -iCj 

i=l,2,3 

Equation  (25)  can  be  transformed  into  natural 
coordinates  by  introducing  the  transformation  matrix  [<I>]  where 
the  coluitvos  of  [<P]  are  the  eigenvectors  found  from  equation 
(31) .  Using  the  transforming  matrix; 


4>2,i 


=  0 


(31) 


’fl 

4'll  4>12  4>i3' 

= 

4'21  4>22  4>23 

^2 

(32) 

"^3 

4»31  4>32  4>33. 

^3, 
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equation  (25)  can  be  decoupled  by  taking  advantage  of  the 
orthogonality  of  the  modal  vectors .  Applying  the 
transformation  matrix  and  premultiplying  by  equation 
(25)  can  be  rewritten  as: 


Assuming  that  the  modal  damping  is  proportional  to  the 
mass  and  stiffness,  equation  (33)  can  be  decoupled  and  written 
as  three  independent  differential  equations  of  motion  as  shown 
below : 


fii+2Co)„iiii+Qfljili=Fi 

(34) 

fi2-^2C0)„2fl2+0)L'n2=-F’2 

(35) 

fi3+2C<.)„3f|3+W^3Tl3=F3 

(36) 

Equations  (34)  through  (36)  were  solved  for 
displacement,  velocity  and  acceleration  using  the  MATLAB™ 
program  provided  in  Appendix  D.  As  in  Chapter  II,  the  critical 
damping  ratio  was  assumed  to  be  0.04.  The  solutions  in  natural 
coordiiiates  was  transformed  back  to  real  coordinates  by 
reapplying  equation  (32). 
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B.  THREE  DOF  FIXTURE  MODEL  RESPONSE 

The  upper  tier  acceleration  response  to  the  half  sine  wave 
impulse  was  modeled  for  the  same  weight  and  stiffness 
configurations  as  presented  for  the  two  DOF  model  in  Figures 
23  through  27  .  The  modeled  acceleration  time  response  and 
shock  spectra  for  five  different  upper  tier  weights  are  shown 
in  Figures  29  through  33  .  As  can  be  seen  from  the  figures,  the 
system  natural  frequencies  for  the  three  DOF  model  are 
approximately  15  to  20  Hz  higher  than  for  the  two  DOF  model 
and  peak  accelerations  are  less.  Figures  34  through  36  show 
the  time  response  and  shock  spectra  comparisons  of  the  two  DOF 
and  three  DOF  models . 

Since  the  anvil  table  is  not  absolutely  unconstrained  due 
to  friction  on  the  bolts  and  contact  with  the  stops,  it  is 
expected  that  the  actual  system  natural  frequencies  will  fall 
somewhere  between  those  of  the  two  DOF  and  three  DOF  models. 
In  the  next  chapter,  the  results  of  the  1/4  scale  model 
testing  will  be  presented  which  will  help  better  predict  the 
response  of  the  actual  fixture. 
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ACCELERATION  (g’s)  ACCELERATION  (g’s) 


40 


MODELED  THREE-DOF  FDCTURE  UPPER  HER  ACCELERATION 


Figure  29 .  Three  DOF  Model  Acceleration  Wave  Fonti  and  Shock 
Spectrum  for  1500  Lb  Upper  Tier  Weight. 
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ACCELERATION  (g’s)  ACCELERATION  (g’s) 


Figure  30.  Three  DOF  Model  Acceleration  Wave  Form  and  Shock 
Spectrum  for  2000  Lb  Upper  Tier  Weight. 
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ACCEl.ERATION  (g’s)  ACCELERATION  (g’s) 


MODELED  THREE-DOF  FDCTURE  UPPER  TIER  ACCELERATION 


5  Foot  Hammer  Drop  (200  g’s) 
Zeta  =  0104 


Upper  Tier  Weigh:  =  2500  lb 
Loivcr-Tier- Weight-  “-2500  Ib- 
Anvil  Weight  =  4500  lb 


TIME  (sec) 


6  C.18 


MODELED  THREE  DOF  FIXTURE  ACCELERATION  SHOCK  SPECTRA 


40  60  SO  1 

.00 

FREQUENCY  (Hz) 

Figure  31.  Three  DOF  Model  Acceleration  Wave  Form  and  Shock 
Spectrum  for  2500  Lb  Upper  Tier  Weight. 
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ACX'EI.ERATION  (g’s)  ACCEl-ERATION  (g’s) 


30 


MODELED  rriREE-DOF  FDCTURE  UPPER  TIER  ACCELERATION 


Spectr\im  for  3  000  Lb  Upper  Tier  Weight. 
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Figure  33.  Three  DOF  Model  Acceleration  Wave  Form  and  Shock 
Spectrum  for  3500  Lb  Upper  Tier  Weight. 
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ACCEl£RAT10N  (g’s)  ACCELERATION  (g’s) 


COMPARING  2  DOF  AND  3  DOF  MODEL  TIME  RESPONSE 


TTME  (See) 


COMPARING  2  DOF  AND  3  DOF  MODEL  SHOCK  SPECTRA 


FREQUENCY  (Hz) 

Figure  34  .  2  DOF  and  3  DOF  Model  Time  Response  and  Shock 


Spectra  Comparisons  for  2500  Lb  Upper  Tier  Weight. 
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ACCELERATION  (g’s)  ACCELERATION  (g’s) 


COMPARING  2  EXDF  AND  3  DOF  MODEL  TIME  RESPONSE 


5  Foot  Hammer  Drop  (200  g’s) 

- 2eta  *  0:04 - 

Upper  Tier  Weight  =  3000  lb 
Lower  Tit-r  Weirjit  =  2700  lb 


\  M 


—  2  DOF 

—  3  DOF 


;  ^  f 


j _ _ ! _ 


:1  L:  1-'  U;  L  -  j 

•  i  ,  J/  I 


0  0.05  0.1  0.15  0.2  a25  0_5  0J5  0.4  0.45 

TIME  (See) 


COMPARING  2  DOF  ANT)  3  DOF  MODEL  SHOCK  SPECTRA 

5  Foot  Hammer  Drop  (200  g’s) 

- 1 - - 1 - Upper-Ticr-Weighv  =-3000-  lb - 

; ;  2  DOF  Lower  Tier  Weight  =  2700  lb 


-.-r * 3  DOf^ 

Ml-. 


FREQLTENCY  (Hz) 


Figure  35.  2  DOF  and  .3  DOF  Model  Time  Response  and  Shock 

SpecLra  Comparison  for  3000  Lb  Upper  Tier  Weight. 
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COMPARING  2  DOF  AND  3  DOF  MODEL  TIME  RESPONSE 

n  1  ■  ^  i  T 

5  Foot  Examiner  Drop  (200  g's)  ,  ' 

. .  Zeta.  =  0.04 .  ,  ]  ; 

Upper  Tier  Weight  =  3500  lb  ' 

Lo'wr  Tier  Weight  2700  lb 
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COMPARING  2  DOF  AND  3  DOF  MODEL  SHOCK  SPECTR.A 


1  ^  ^  5  Foot  Hammer  Drop  (200  g's) 

I  "  _ Upper  Tier  Weight  «  3500  lb 

!  Lower  Tier  Weight  >“  2700  lb 

i'ODOF'"  “  i . . 
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FigT.ire  36.  2  DOF  and  3  DOF  Model  Time  ReGpo.nse  and  Shock. 

Spectra  Cerpar i.son  for  3500  LL  Upper  Tier  VJeight  . 
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V.  1/4  SCALE  MODEL  TESTING  AND  RESULTS 


To  verify  the  mathematical  simulation,  a  1/4  scale  model 
of  the  fixture  was  fabricated  for  shock  testing.  The  principal 
intent  of  the  testing  was  to  confirm  the  analytical  methods 
employed  were  appropriate  for  the  application  and  that  the 
expected  frequencies  characteristics  could  be  obtained.  The 
following  sections  describe  the  test  procedure  used  and  the 
results  obtained. 

A.  EXPERIMENTAL  SET  UP 

To  supply  the  inpulse  excitation  to  the  model,  a  drop 
table  was  constructed  which  allowed  dropping  the  fixture  model 
onto  a  load  cell  from  various  heights.  The  range  of  the  load 
cell  was  0-5000  lbs  with  a  sensitivity  of  0.106  mv/lb.  The 
upper  tier  accelerometer  had  a  range  of  0-500  g's  with  a 
sensitivity  of  9.98  mv/g.  Using  a  Hewlat  Packard  3562A  D^Tiamic 
Signal  Analyzer,  the  transfer  function  between  the  input 
excitation  and  upper  tier  acceleration  was  obtained.  The 
experimental  set  up  is  shown  in  figure  37  and  a  schematic  of 
the  instrumentation  is  shown  in  figure  38. 

B.  1/4  SCALE  MODEL 

The  fixture  model  was  built  as  clo.sely  as  possible  to  a 
1/4  scale  version  of  the  proposed  fixture.  Ebcceptions 


included  using  aluminum  instead  of  steel  for  the  spring  I- 
beams  and  an  additional  1/4  scale  reduction  in  the  weight  of 
the  tiers.  The  upper  tier  weight  as  tested  was  43.5  lbs,  the 
lower  tier  weight  was  55.2  lbs  and  the  anvil  table  was  109.1 
lbs.  The  per-beam-stif 'ness  of  the  aluminum  I-beams  was 
2.09x10*  Ibs/in. 

To  help  verify  the  boundary  conditions  used  in  the 
mathematical  model,  the  spring  beams  were  mounted  using  the 
scime  clamping  design  proposed  for  the  full  scale  prototype. 
The  following  exaimple  provides  the  dimensions  and  calculations 
for  determining  the  expected  system  natural  frequencies  for 
the  three  DOF  model . 


Example:  Detemiine  System  Natural  Frequencies  For  1/4  Scale  Model  Fixture: 
A.  SPRING  STIFFNESS: 

The  spring  stiffness  is  based  on  the  following  beam  model; 

y 

A 

Simple  Supports  -  Symetric  Loads 


1)  Dimensions: 

Length: 

Load  Applicaton  Point; 
Area  Moment  of  Inertia: 
Young's  Modulus: 


Upper  Tier 

L  ]  =12.5  in 
a  I  =  3.25  ui 
I  I  =0.045  in* 
E  =  10  10*  psi 


Lower  Tier 

L  2  =12.5  in 
aT  =3.25  in 
1  2  -  0  045  in* 
E  =10  10*  ps; 
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2)  Stiffness: 


From  the  beam  bending  equation: 


Let  x=a  and  F=W/2; 


W'fl 

yAB* - ^  (4  a-  3  L) 

^  12EI 


Solving  for  the  stiffness: 


,  12  E  I  ]b 

yAB  a^(3L-4a)  ^ 


For: 


Upper  Tier 


Lower  Tier 


(Convert  inches  to  feet)  KB  j 


12E  I 


a  j  i  3  ■  L  I  -  4  a  j  i 


12  KB 


12  E  K 


a'>  ;3L2“4'a'^'. 


(per  beam) 

K3  ■  =2.504-10*- 
tt 

KBt  =2.504- 10^- 
ft 

Upper  Tier 

Lower  Tier 

Input  Number  of  Beams  Per  Tier:  Beams  j  =  2 

Beams 2 

K  j  =  KB  j  Beams  • 

Kt  =  kb  2  Beams-) 

Tier  Stiffness:  K,=5.008-lo’  - 

'  ft 

K,  =5.008-10’  - 

ft 

B.  SYSTEM  MASSES;  (Includes  Effective  Sprina  Mass) 

Upper  Tier  Weiaht:  (W'^) 

Intermediate  Weiaht:  (Wj) 

Anvil  Table  Weiaht:  IV7-^) 

W,  =43.5 

W2  =55.2 

W3  =  109  1 

M  ,  = — !- 
‘  32.3 

W2 
= — =- 

^  32.2 

W3 

M  1  =  — 

32.2 

M,  =1.351  lb— - 
‘  .2 

M2  ==1714  lb -5-- 

M3  =  3  388  lb 

sec 


sec 


ft 

sec‘ 
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C.  TIER  NATURAL  FREQUENCIES: 


D.  System  Natural  Frequencies: 

Defining  B  as; 


Upper  Tier 


F  =J_. 

F  ,  »  96.904  Hz 


Lower  Tier 


“  2rzjM2 


F  2  =  86.023  Hz 


B  =■ 


K,  K,  K , -K, 


M  1  M3 


M., 


,  i|  1  I  ,  4.K,K,  I 

f  .  =J_.  il|B-  !B‘ - 1 — i_.;m  , I 

ni  I  1  V/  v.<  v.<  \  1  -  ■’/  i 


First  Frequency; 


2  z  12! 

‘  M  ; 


f^l  =77.762  Hz 


N 


M  [  M  2  M  3 


'  02 


I  ji  u  i'p2  ^•^r^2 


IB-  jB‘- 

J  M1M2M3 


.  (M  1  -  M  2  “  M  3^,  i 


Second  Frequency:  fj^  =  14^  y44  Hz 


C .  EXPERIMENTAL  RESULTS 

The  first  experiment  performed  was  determining  the  static 
stiffness  of  an  assembled  tier  containing  three  spring  beam 
elements.  Using  the  simply  supported  beam  model,  the  expected 
stiffness  of  the  tier  was  calculated.  The  tier  was  then  place 
on  a  load  machine  and  loaded  up  to  lOOC  lbs  while  measuring 
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the  static  deflection.  Figure  39  shows  the  analytical 
stiffness  versus  the  experimental.  From  the  plot,  the  actual 
stiffness  was  measured  to  be  43,500  Ib/in  which  is  about  5 
percent  higher  than  the  analytical  stiffness.  This  result 
tends  to  justify  the  simply  supported  boundary  conditions  used 
in  the  numerical  model. 


EXPERIMENTAL  VS.  CALCULATED  TIER  STIFFNESS 

0.02S 

0.02 

0.013 

Displacement 
(Inches) 

0.01 

0.005 

0 

200  300  400  500  600  700  *00  900  lOOO  1100 

LOAD  (Lbs) 

Figxtre  39.  Experimental  Tier  Stiffness  versus  Analytical. 

The  model  was  then  installed  on  the  drop  table  to 
determine  its  response  to  a  six  inch  table  drop.  Using  the 
load  cell  as  the  input  and  upper  tier  acceleration  as  the 
output,  the  tremsfer  function  for  a  ten  drop  average  was 
computed  using  the  HP  3562A.  The  acceleration  data  was 
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filtered  using  an  exponential  window  with  a  period  T=400  msec. 
The  data  collection  was  triggered  off  the  load  cell  input 
using  a  -5  msec  trigger  delay.  The  transfer  function  with 
phase  shift  is  shown  in  Figure  40.  The  experimental  transfer 
function  reveals  distinct  peaks  at  80  Hz  and  160  Hz  with 
associated  180  degree  phase  shifts.  These  results  agree  very 
closely  with  the  mathematical  predictions  of  fni  =  77.8  Hz  and 
f„2  =  148  Hz. 

Figure  41  shows  the  \infiltered  load  cell  and  acceleracion 
output  sampled  at  10,000  Hz.  The  load  cell  data  shows  a 
maximum  impulse  of  approximately  35,000  lbs  with  a  duration  of 
about  10  msec.  The  corresponding  upper  tier  acceleration  is 
17  0  g's.  Usi.ng  a  logarithmic  decrement  as  an  approximation, 
the  critical  damping  ratio  was  determined  to  be  two  percent. 

The  experimental  results,  although  not  extensive,  support 
the  mathematical  model  very  well.  They  demonstrate  that  using 
the  proposed  fixture  will  transform  the  high  frequency,  high 
impact  energy  of  the  hcimmer-anvil  impact  into  discrete 
equipment  excitation  frequencies,  tunable  to  shipboard 
conditions.  With  full  see Ve  fixture  testii  g  beginning  in  July 
1993,  promising  results  are  exjjected. 
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ACCELKRATION  (g’s)  FORCE  (Lbs 


I 


Accelerometer  Output.  Sample  Frequency  10,000  Hz. 
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VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  mathematical  model  and  experimental  results  provided 
in  this  study  demonstrate  great  promise  in  the  concept  of 
using  a  tuned  two  DOF  equipment  mounting  fixture  on  the  MWSM. 
Using  frequencies  provided  by  NSWC,  the  mathematical  model 
provided  the  required  stiffness  and  mass  parameters  necessary 
to  design  the  fixture.  The  design  was  then  scaled  to  one 
quarter  and  an  experimental  model  constructed  for  impact 
testing.  The  results  of  the  model  testing  agreed  extremely 
well  with  the  results  predicted  from  the  mathematical 
modeling.  The  following  is  a  list  of  specific  comments 
concerning  the  tuned  fixture: 

-  The  fixture  designed  in  this  study  is  not  expected  to  be 
an  all  inclusive  tuned  fixture  capable  of  testing  all 
mediumweight  equipment.  Its  specific  design  is  to  provide 
excitation  at  30  Hz  and  80  Hz  for  equipment  up  to  3500 
lbs.  As  more  information  becomes  available  on  localized 
ship  structure  excitation,  a  series  of  tuned  fixtures  can 
be  constructed  (economically) ,  each  designed  for  a 
specific  location  on  the  ship. 

-  Since  the  weight  capacity  of  the  MT-JSM  is  limited  to  7  4  00 
lbs,  equipment  weight  using  a  two  DOF  fixture  will  be 
limited  to  approximately  3500  lbs.  The  ratio  between  the 
upper  and  lower  tier  mass  is  limited  requiring  them  to  be 
nearly  equivalent  in  magnitude.  From  the  DTRC/UERD 
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predictions,  this  should  not  be  a  significant  problem 
since  heavier  equipment  will  not  generally  require  a  two 
DOF  fixture  to  simulate  shipboard  shock. 

-  The  I-beam  spring  design  will  work  very  well  for 
frequencies  above  25  Hz,  but  would  not  be  suitable  for 
lower  excitation  frequencies.  Since  the  spring  beams  are 
limited  to  approximately  sixty  inches  in  length  on  the 
MVJSM,  the  vertical  tier  displacements  required  for  lower 
frequencies  would  result  in  excessive  bending  stress  in 
the  I-beams.  For  this  reason,  it  is  recommended  that  the 
coil  spring  design  be  reinvestigated  tor  lower 
frequencies . 

-  On  completion  of  the  prototype  fixture  in  June  1993,  a 
series  of  calibration  tests  will  be  performed  to  determine 
the  actual  response  of  the  fixture.  Based  on  these  test, 
the  mathematical  model  can  be  updated  with  structural 
dcimping  characteristics  as  well  as  adjusting  the  boundary 
conditions  if  necessary.  This  new  model  will  then  be 
available  for  the  design  of  future  fixtures. 

The  use  of  the  two  DOF  equipment  mounting  fixture  for 
shock  qualifying  equipment  on  the  MWSM  is  in  line  with  the 
Navy's  goal  of  shock  hardening  ships.  Incorporating  this 
concept  into  the  current  shock  qualification  procedures  w’ill 
provide  increased  reliability  of  essential  equipment  subjected 
to  shock  from  underwater  explosions. 
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APPENDIX  A.  TWO  DOF  TUNED  FIXTURE  PROGRAM 


%  Program  TWODOF.M 

%  Lt .  David  M.  Cox 

%  Naval  Postgraduate  School 

%  This  MATLAB  program  plots  the  relative  displacement, 
%relative  velocity  and  absolute  acceleration  of  a  TDOF  mass, 
%spring,  damper  system  subjected  to  base  excitation  (half-sine 
%pulse) .  Additionally,  the  shock  spectra  of  an  undamped  SDOF 
%system  is  calculated  using  the  TDOF  upper  tier  acceleration 
%as  the  base  excitation  to  the  SDOF  system. 

clear 


%User  Defined  Variables 


G=203 

Wtl=3500 

Wt2=2700 

kb=2 . 6e6 

kl=3‘kb; 

k2=4*kb; 

2eta= . 04 


%  Acceleration  Magnitude  in  g ' s 
%  Weight  of  Upper  Tier  in  Lbs 
%  Weight  of  Lower  Tier  in  Lbs 
%  Stiffness  per  Beam 
%  Upper  Tier  Stiffness  in  Lbs ft 
%  Lower  Tier  Stiffness  in  Lbs/ft 
%  Critical  Damping  coefficient 


%Formu].ate  the  Base 
dt=0 .0001 
w=1000*pi 
N=5000; 

T=N*dt 

t=[0:dt: (N-1) *dt] 
for  i=l:N 
if  t(i)<=0.001 

anvil(i)  =  G*32 .2*sin(w*t 
else 

anvil (i)  =  0  ; 

end 

end 


Acceleration 

;  %  time  step  for  sampling 
;  %  base  excitation  frequency 
;  %  Number  of  steps 
;  %  Total  time  record 
;  %  Time  vector  from  t=0  to  t=T 


%  Loop  to  produce  half -sine  wave  pulse 


i )  ) 


%Calculate  the  mass  and  mass  ratio 
ml=Wtl/32.2  ; 
m2=Wt2/32.2  ; 
m= [ml , 0 ; 0 , m2 ] ; 
alpha=ml /m2 ; 
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%Calculate  Circular  Tier  Natural  Frequencies 
wl=sqrt (kl/ml) ; 
v2-so''t  (k2/m2 )  ; 

%Systein  Natural  Frequencies 

var=sqrt  (  (v;1^2  +  alpha*wl^2  +  w2''2)^2  -  (2*wl*w2  )  ^^2 ) 
wnl=  ( 1/sqrt  (2  )  )  *sqrt  (wl"^2  +  alpha*wl^2  +  w2^2  -  var) 
wn2=  ( 1/sqrt  ( 2  )  )  *sqrt  (wl''2  alpha*wl^2  +  w2''2  +  var) 
%Uncouple  the  equations 
ull=l; 

u21=(kl  -  wnl''2  *  ml)/kl; 
ul2=l; 

u22=(kl  -  wn2^2  *  ml)/kl; 

U-[ull,ul2;u21,u22]  ; 

%Decouple  the  system 
M=U' *m*U 
%£tiffness  Matrix 

k= [kl, -kl;-kl,kl+k2] ;  %  coupled  stiffness 

K=U' *k*U 

%Force  Coefficients 

FC=(U' * [-ml; -m2] ) 

%Solve  the  uncoupled  Equations  Of  Motion 

%  {nddotl }  +  [2  *zeta*Vi'nl  ]  *  {ndotl }  +  [wnl^2  ]  *  {nl }  =FC11* anvil  /Ml 
%  {nddot2 )  +  [2  *zeta*wn2 ] * {ndotl} + [ wn2  ^2 ] * {n2 } =FC21*anvil /M2 

%State  Space  Matrix  for  Equation  #1 

A1=[0, 1;  -  (wnl'^2)  ,  -2*zeta*wnl] 

B1=[0;FC(1;1) /M(l:l) ] ; 

Cl=tl,0;0,l;-  (wnl''2  )  ,  -2*zeta*vml]  ; 
Dl=t0;0;FC(l:l)/M(l:l) ) ; 

[yl] =lsim(Al,Bl,Cl,Dl, anvil, t) ; 

%State  Space  Matrix  for  Equation  #2 

A2  =  [  0 , 1 ;  -  {wn2''2  )  ,  -2*zeta'*wn2  ] 
B2=[0;FC(2:2)/M(4:4)]; 

C2  =  [  1 , 0  ;  0 , 1 ;  -  (wn2''2 )  ,  -  2*zeta*vm2  ]  ; 

D2=[0;0;FC(2 :2) /M(4 :4) ] ; 

[y2 ] =lsim(A2 , B2 , C2 , D2 , anvil , t ) ; 

%Couple  The  Equations 

Y1DIS= (ull*yl { ; , 1)  +  ul2*y2( ; , 1) )*12; 

YlVEL=ull*yl ( : ,2)  +  ul2*y2{:,2); 

YlACC=ull*yl( :,3)  +  ul2 *y2 { ; , 3 ) ; 

%Absolute  acceleration  in  g's 

X1ACC= (YlACC+anvil' ) /32.2; 

%Convert  to  in. 

Y2DIS= {u21*yl ( : , 1)  +  u22*y2(:,l))*12; 

Y2VEL=u21*yl ( : , 2)  +  u22*y2(:,2); 

Y2ACC=u21*yl { : , 3)  +  u22 *y2 ( : , 3 ) ; 

X2ACC= { Y2ACC+anvil ’ ) /32 . 2 ; 
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plot ( t, XIACC) , grid 

title ( 'MODELED  TWO  DOF  FIXTURE  UPPER  TIER  ACCELERATION') 

xlabeK 'TIME  (sec) ' ) 

y label ( 'ACCELERATION  (g''s)') 

gtext ( ' 5  Foot  Hammer  Drop  (200  g''s)') 

gtext( 'Upper  Tier  Weight  =  3500  lb') 

gtext( 'Lower  Tier  Weight  =  2700  lb') 

gtext ('Zeta  =  0.04') 

meta  fig28.met 

plot (t, YIDIS, t, Y2DIS) .grid 
xlabel ( ' TIME { Seconds ) ' ) 

ylabel  ( 'RELATIVE  DISPLACEMFAT'  (Inches)') 

title ( 'MODELED  TIER  DISPLACEMENT  RELATIVE  TO  THE  ANVIL  TABLE' ) 

gtext (' Hammer  Height  =  5  feet  (200  g''s)') 

gtext ('1  -  Upper  Tier  Weight  =  3500  lb') 

gtext ('2  -  Lower  Tier  Weight  =  2700  lb') 

gtext('Zeta  =  0.04') 

gtext ( ' 1 ' ) , gtext ( ' 2 ' ) 

meta  fig20.met 

pause 


%Calculate  the  shock  spectr 
NF=120; 

DF=:1; 

SF=1  ; 

Freq=(SF;DF:NF] ; 
wn=2*pi*Preq; 

F=[0;-1] ; 

GG=tl, 0] ; 

H=l-1]  ; 


%  Number  of  Frequencies 
%  Frequency  Increment 
%  Start  Frequency 
%  Frequency  Vector 

%  State  Space  Matrices 


for  i=l;NF 


%  Step  Thru  Natural  Frequencies 


E=[0,l;-(wn(i)'2),0]; 

GG=[wn(i)  ''2,0]; 

[yspec] =lsim(E, F, GG, H, XIACC, t ) ; 
xspec  =  yspec  +  XIACC; 
maxspec ( i ) =max ( abs ( xspec ) )  ; 
end 


plot (Freq, maxspec ) 

title ( 'MODELED  TWO  DOF  FIXTURE  ACCELERATION  SHOCK  SPECTRA') 
xlabel ( ' FREQUENCY  (Hz) ' ) 
ylabel ( 'ACCELERATION  (g''s)') 
grid 

gtext ('fnl  =  38  Hz') 
gtext ('fn2  =  78  Hz') 
meta  fig28.met 
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APPENDIX  B.  TUNED  FIXTURE  DESIGN  DRAWINGS 


DESIGN  DRAWINGS 


ITEMIZED  PART  LIST 


1  Piece  Number 

Item 

Size  (inches) 

No.  Required 

'  1 

Shipbuilding  Channel 

Ixll.lU 

2 

'l  2 

I-Beam 

W5  X  18.5# 

4 

i  3 

I-Beam 

W8  X  67# 

3 

1  ^ 

Block 

l"x2  1/4"  X  10" 

28 

1  5 

Stiffner 

See  Figure  2  Sheet  2 

6 

6 

I-Beam 

W4x  13# 

7 

1  ^  - 

Stiffher 

See  Figure  2  Sheet  5 

12 

1  8 

Block 

I"x2"x5  1/8" 

28 

1  9 

Clamp 

See  Figure  4  Sheet  1 

14 

1  10 

Clamp 

See  Figure  5  Sheet  1 

28 

i _ y _ 

Washer 

1  l/8"IDx2"OD 

84 

^ m. 

Hexagon  Head  Bolt 

l"-8x4  1/4"  Long 

84 

Notes : 

1)  Top  flange  of  piece  number  1  shall  be  burned  off  or 
cut  off  to  a  width  of  1  3/4". 

2)  Piece  numbers  9  and  10  will  be  shaped  to  fit  inner 
surface  of  piece  nunbers  1  and  6  respectively. 

0)  Hexagon  head  bolts  may  be  replaced  with  eq\iivalent 
alien  head  bolts  ease  installation  of  piece  10  to 
piece  2 .  Depending  on  weld  thickness  and  assembly 
tolerance,  use  of  a  socket  to  tighten  hexagon  head 
bolts  may  be  difficult. 
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APPENDIX  C.  MATHCAD®  WORKSHEETS 


APPENDIX  C1 

MathCad  Worksheet  To  Find  Two  DOF  Fixture  Ucoupled  Frequencies: 


input  the  Desired  Coupled  Natural  Frequencies:  (Hz) 

=30.00  -- 80.00 

Input  The  Mass  Ratio: 

Mass  Ratio:  a  =  I 

Initial  Guess  at  uncoupled  frequencies:  f  ]  =40  f  2  =40 

Given  the  initial  guess  for  and  f2,  the  following  equations  are  iteratively  solved  for  the  required 
values  of  fi  and  f2. 

Given 


SOLUTION 

Required  Uncoupled  Tier  Natural  Frequencies:  f  ]  =  49. %3  f  2  =  48  036  Hz 
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APPENDIX  C2 


Mathcad  Worksheet  to  Determine  System  Natural  Frequencies  Of  Two  DOF  Fixture: 


A.  SPRING  STIFFNESS: 


The  spring  stiffness  is  based  on  the  following  beam  model: 

y 


Supply  the  following  Dimensions; 


1)  Dimensions: 

UoDerTier 

Lower  Tier 

Length: 

L  I  =  so  in 

L  2  =  50  in 

Load  Application  Point: 

a  I  =  13  in 

8  2=13  in 

Area  Moment  of  Inertia; 

1,  ^113  ui* 

I2  =113 

Young’s  Modulus: 

E  =30  10*  psi 

E  =30  10*  psi 

2)  Stiffness: 


From  the  beam  bending  equation;  v  (x^  +  3  -  3  l) 

'  ^  6  El 


Let  x=a  and  F=W/2; 


W 

>AB= — ^•(4a-3L) 
^  12E1 


Solving  for  the  stiffness: 


W  _  12  E1  lb 

>AB  B^(3L-4a) 


For; 

UoDer  Tier 

12  El  , 

(Convert  inches  to  feet) 

KB,  ^ - 12 

a  ,  (3  1,  -  4  a  j) 

Lower  Tier 

12  El  , 

KB2  = - = -  12 

82  (^'1^2” 


(per  beam) 


KB  1  =2.947-10*  - 
‘  ft 


KB  9  =2.947-10* 

ft 
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Input  Number  of  Beams  Per  Tier; 


Beams  |  =  2 


Beams  2  =  ^ 


K  j  =  KB  ]  Beams  j 


Tier  Stiffness:  K, -5  895-10*  - 

*  ft 


K  2  =  KB  2  Beams  2 

K,  =  1.179-10'^  - 

^  ft 


B.  SYSTEM  MASSES:  (Includes  Effective  Sjirina  Mass 
Upper  Tier  Wieght:  W  j  =  2500  lbs  Lower  Tier  Weight:  W  2  =  2500  lbs 


*  IT 


M  I  =77.64  lb 


M-,  =77.64  lb 


Mass  Ratio  a: 


C.  TIER  NATURAL  FREQUENCIES: 


.ower  Tier 


f, 


1  K, 


2-rc  JM 


D.  SYSTEM  NATURAL  FREQUENCIES: 


f  -  J_.  ^ 

^  2n^M2 


fj  =43.855  Hz  f2  =62.02  Hz 


fnl  (f2)^j  -  (21,  (2)' 

V2 


f„,  =  33.565  Hz 


fn2  ^■7’-](^l)^+ “  (^2)^+^*  (f  1.'^^- “  v^ l)^+ (^2)^1  -{^'^vhy 

V2 


fjj2  =81.033  Hz 


APPENDIX  C3 


Mathcad  Worksheet  to  Find  Two  DOF  Coupled  System  Natural 
Frequencies  as  a  Function  of  the  Uncoupled  Tier  Frequencies 


Pick  a  Value  For  the  Lower  Tier  Uncoupled  Natural  Frequency  (F2);  F2  -  48  Hz 
Specify  the  Mass  Ratio;  a  =  1  0 


Solve  for  the  Coupled  System  Natural  Frequencies  as  a  Function  of  F1: 

k  ^  I  .  100  Flj^  =k 


FNl. 


r 


/t 


iiFl.i 

V  ,  k- 


a 


■FI  ‘ 
I  k 


(F:r 


iFl 


s2 


a  FI. 


(F2)‘ 


FI  •F2'',^ 


FN2j^  =■ 


liFL  , 
'V  '  k, 


\2 


a  (Fl^)^+  (F2r- j[(Fg^+a  ;F1^^-(F2)^ 


-i2Fl,F2V 


i  =40.60 


PICK 

UPPER  TIER  COUPLED 


FI 

t 

12 

41 


12 

44 

45 

46 

47 

48 

12 

E 

51 

52 

53 

54 

22 

56 

57 

22 

59 


FNl  FN2. 

_ I  _ i 


27.936 

68.728 

28  198 

69.792 

28.446 

70.872 

28.679 

7 1 .969 

28.899 

73.081 

29  108 

74.207 

29.304 

75.348 

2949 

76.501 

29.666 

77.666 

29.832 

78.842 

29.989 

80.029 

30.138 

81.226 

30.279 

82.433 

30.413 

83.648 

30.54 

84.872 

30.661 

86.104 

30.775 

87.3'44 

30884 

88  59 

30.987 

89.843 

31.086 

91.103 

31.18 

92.368 

COUPLED  SYSTEM  NATURAL  FREQUENCIES 

a  =  1  F2  =  48 


UPPER  TIER  NATURAL  FREQUENCY 
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APPENDIX  C4 


MathCad  Worksheet  To  Find  Three  DOF  Fixture  Ucoupled  Frequencies: 


Mass  Ratio: 


Input  the  Desired  Coupled  Natural  Frequencies:  (Hz) 

f„,  ^30.00  -80  00 

Input  The  System  Tier  Weights:  (lbs) 


1  =2500  lbs 

W2  =2500  lbs 

W3  =4500 

Wi 

Wo 

W  , 

1 

P  ^  — 

Y  = - 

W3 

a  =  1 

3  ^0  556 

Y  =  0  556 

Supply  the  Program  With  An  Initial  Guess  at  f-i  and  (3: 


f;  =40  Hz  f2  =40  Hz 


Given  the  initial  guess  for  fi  and  (2,  the  following  equations  are  iterativley  solved  for  the  required 
values  of  f-|  and  f2 


l  \ 


ff,)^-(l  +  a)+ (f2)'(U(3) 


a)+,f2;  (1  +  P)  -4 


Y+  1) 


\f2 


=  Findff  1^2! 


Required  Uncoupled  Tier  Natural  Frequencies:  f  j  =  54  077  Hz 


f2  =  30.545  Hz 
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APPENDIX  C5 


Mathcad  Worksheet  to  Determine  System  Natural  Frequencies  Of  Three  DOF  Fixture: 


A.  SPRING  STIFFNESS: 

The  spring  stiffness  is  based  on  the  following  beam  model; 

y 


Simple  Supports 

'  Symetric  Loads 

Supply  the  following  Dimensions: 

1)  Dimensions: 

Upi 

oer  Tier 

Lower  Tier 

Length: 

^1 

=  50 

in 

L2  ^50 

in 

Load  Application  Point: 

=  13 

in 

&2  ^'3 

in 

Area  Moment  of  Inertia: 

h 

=  11.3 

in 

I2  =113 

4 

in 

Young's  Modulus: 

E  : 

=  30  10* 

psi 

E  =30  10’ 

*  psi 

2)  Stiffness: 

From  the  beam  bending 

equation: 

>  AB 

6  El  ' 

x^+  3  a^-  3  l) 

_  W 

(4  a-  3  L) 

Let  x=a  and  F=W/2; 

>■  AB 

12E  1 

..  W 

12E  1 

lb 

Solving  for  the  stiffness: 

K=- 

y 

AB 

(3  L-  4  a) 

in 

(Convert  inches  to  feet) 


Upper  Tier 


a  1^  (3  L  1  -  4  a  ,  ) 


KB,  :: - 


Lower  Tier 

12EI, 


^  ®2' 


(per  beam) 


KB  ,  =2  947-10*- 
‘  ft 


KB  2  =  2  947-10* 


I 


Upper  Tier 

Input  Number  of  Beams  Per  Tier:  Beams  j  =2 


Tier  Stiffness: 


K  ]  =KB  )  Beams  j 
K  ,  =5.895- iO'^  - 

'  n 


B.  SYSTEM  MASSES:  (Includes  Effective  Spring  Mass) 
Upper  Tier  Weight:  (W-j )  Intermediate  Weight: 


^'1 

=  2500 

W.  =2500 

w. 

M  , 

M.  = — ^ 

1 

32.2 

--  32  2 

M  , 

=  77.64  Ib-iL 

M  .  =77.64 

sec 


sec 


C.  TIER  NATURAL  FREQUENCIES:  Upper  Tier 


IK 


f 


1  •; 


2  Tt  J  M  j 

f,  =43.855  Hz 


D.  System  Natural  Frequencies: 

Defining  B  as. 


Lower  Tier 


Beams  2  - ^ 


K2  =  KB  2  Beams  2 

K-,  =  1179-10'^  - 

ft 

Anvil  Table  Weight:  (W->) 

W3  =4500 

W-, 


M 


32.2 


M3  =  139.752  lb 


Lower  Tier 


ft 

scc‘ 


IK. 


2  n 


f2  =62.02  Hz 


K  ,  Kj  K  1  +K. 
B  =— + 


M I  M3  M2 


f. 


f 


:-L.  l.|B- 


2«^|2[  ^ 

First  Frequency:  fjjl=44  64  Hz 


2  4K,K2  ,1 

^  M  M  (Mi^M2^M3)I 

M  ]  M  2  M  y  ^  I 


*n2  1:^ 

2n  : . 


B-t- 


B- 


4.K,K2  , 


M  ]  M  2  M  3 


^M  ]  +  M  2  -t-  M  3  ': 


Second  Frequency:  1^=88  528  Hz 
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APPENDIX  C6 


Mathcad  Worksheet  to  Find  Three  DOF  Coupled  System  Natural 
Frequencies  as  a  Function  ot  the  Uncoupled  Tier  Frequencies 

Pick  a  Value  For  the  Lower  Tier  Uncoupled  Natural  Frequency  (f2)  f  o  =  48  Hz 


Input  The  System  Tier  Wieghts 

W  1  =  2500  lb 

V'2  =2500  lb 

V,':,  =4500 

and  calculate  mass  ratios; 

n  '^2 

V-’l 

a  =  — 

P  r - 

Y  " - 

W2 

^•'3 

^'3 

a  =  1 

P=0  556 

Y  =  0  556 

Solve  for  the  Coupled  System  Natural  Frequencies  as  a  Function  of  ri: 

Define  Mass/Stiffness  Constant:  B.  -,fi  ( I  +  a)  *  (f '>’1^  ( I  +  p) 

‘I  k, 


-k 


-  4 


I'f  I  1 

\  k.' 


(pTy+  1) 


I 


3  il  |r  ,  4./f,  ^-.(f2i’(P+Y*  l)j 

-J2  [  ^  ^  \  ‘k.'  ^  • 


P/CK  i  -  30..  50 

UPPER  TiER  COUPLED 


n.  ^nl.  ^ti2. 

I  I  I 


30 

31.593 

66.226 

31 

32.383 

66763 

32 

33  147 

67  33 

33 

33.882 

67.926 

34 

34.59 

68  553 

35 

35  269 

69.21 1 

36 

35.92 

69  898 

37 

36.542 

70  616 

38 

37  136 

71  365 

39 

37.703 

72  142 

40 

38.242 

72949 

38  75J 

n  784 

42 

39.24 

74647 

43 

39.702 

75.537 

44 

40  138 

76452 

45 

40.552 

77.392 

46 

40.943 

78.356 

47 

■3’  313 

79.342 

48 

yy.l 

30.35  r 

49 

■]  ^^2 

SI. 379 

50 

|4"  .06 

32.427 

COUPLED  SYSTE?ifl  NATURAL  FREQUENCIES 

a  =  l  f2=48 


150, 

y 

FN? 

^  ■ 

y 

— 

FN1 

_ Q 

0  20  40  60  *0  100 


f>k 

UPPER  TIER  NATURAL  FREQUENCY 
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APPENDIX  D.  THREE  DOF  TUNED  FIXTURE  PROGRAM 


%  Program  THP.ELDOF.M 
%  Lt .  David  M.  Cox 
%  Naval  Postgraduate  Sch  >ol 


%  This  program  plots  the  relative  displacement , relative 
%velocity  and  absolute  acceleration  of  a  three  DOF  mass, 
%spring,  damper  system  subjected  to  base  excitation  (half-sine 
%pulse) .  Additionally,  the  shock  spectra  of  an  undamped  SDOF 
%system  is  calculated  using  the  upper  tier  acceleration  as  a 
%base  excitation  to  the  SDOF  system. 


%User  Defined  Variables 


HH=5  % 
Wtl=3000  % 
Wt2=2500  % 
Wt3=4500  % 
kbean=2 . 947e6  % 
kl=3*kbeam  % 
k2=4*kbeam  % 


2eta= . 04 


g=32.2 


% 


Drop  height  in  feet 
VJeight  of  Upper  Tier  in  Lbs 
Weight  of  Lower  Tier  in  Lbs 
Weight  of  Anvil  Table  in  Lbs 
Stiffness  of  one  Beam  Lbs/ft 
Upper  Tier  Stiffness  in  Lbs/ft 
Lower  Tier  Stiffness  in  Lbs/ft 
Critical  Damping  coefficient 
acceleration  of  gravity 


%Fonnulate  the  Base  Acceleration 


%  time  step  for  sampling 
%  Number  of  steps 
%  Velocity  of  Table  after  impact 
%  From  Figure  19 
%  Total  time  record 
%  Time  vector  from  t=0  to  t=T 


dt=0.-0001; 

N=5000; 

Tblvel=sqrt ( 2 *g*HH) ; 

PkAccel  =  ll . 3  5*Tblvel 
T=N*dt; 

t=[0:dt: (N-1) *dt] ; 
w=1000*pi ; 

for  i=l:N  %  loop  to  produce  half-sine  wave  pulse 

if  t(:.)<  =  0.001 

Force(i)  -  (Wt3 ) *PkAccel*sin (w* t { i ) ) ; 
else 

Force (i)=  0; 


end 

end 


%Calculate  the  mass  and  mass  ratio 
ml=Wtl/32.2  ; 
m2=Wt2/32.2  ; 
m3=Wt3/32.2  ; 
m=  [  ml ,  0 , 0  ;  0  ,  in2 , 0  ;  0 , 0  ,  m3  ]  ; 
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alpha=ml/m2 ; 
beta=m2  /it3  ; 
gainma=ml/m3  ; 

%Calculate  Circular  Tier  Natural  Frequencies 
wl=sqrt(kl/rT\l); 
w2=sqrt  (k2/in2)  ; 

%Systein  Natural  Frequencies 

b=kl/iul  +  k2/m3  +  (kl  +  k2 ) /iti2  ; 
wnl  =  0 

wn2=sqrt  (0 . 5*  {b-sqrt  (b''2-  (4*kl*k2)  *  (ml+m2+m3  )  /  (inl*m2*ni3  )  )  )  ) 
wn3=sqrt  (0.5*  (b+sqrt  (b^2-  (4*kl*k2 )  *  (ml+in2+m3  )  /  (inl*in2 *m3  )  )  )  ) 

%Uncouple  The  Equations 
ull=l . 0; 
u21=1.0; 
u3l=l . 0; 
ul2=1.0; 

u22  =  ( kl-ml*vm2  ^2  )  /  kl  ; 

u32=(-kl  +  (kl  +  k2  -  m2  *  vm2^2 ) *u22 ) /k2 ; 
ul3=1.0; 

u23= ( kl-ml*wn3  ^2 ) /kl ; 

u33=(-kl  +  (kl  +  k2  -  m2  *  wn3 '‘2 )  *u23  ) /k2  ; 

U= [ull , ul2 , ul3 ;u21 , u22 , u23 ; u31 , u32 , u33 ] 

%Decouple  the  System 
M=U' *m*U 

%Stiffness  Matrix 

k= (kl , -kl , 0 ; -kl , kl+k2 , -k2 ; 0, -k2 , k2 ] ; 

K-U' *k*U 

%Force  Coefficients 
FC=U ' * ( 0 ; 0 ; 1 ] 

%Solve  the  Uncoupled  Equations  of  Motion 

%  {nddotl)  +  [2*zeta*wnl"2]  *{ndotl)  +  [wnl''2]  *{nl}=FCll*Force/Ml 
%  {nddot2  }  +  (2*zeta*wn2^2  )  *  {ndot2  }  +  [wn2''2  ]  *  {n2  }  =FC21*Force/M2 
%  {nddot3)+ [2*zeta*wn3^2] *{ndot3)+ [wn3^2] * {n3 ) =FC31* Force /M3 

%State  Space  Matrix  for  Equation  #1 
Al= [0 , 1 ; - (wnl^2 ) , -2*zeta*wnl] ; 

B1=[0;FC(1:1> /M(l:l) ] ; 

Cl= [ 1 , 0 ; 0 , 1 ; - (wnl^2 ) , -2*zeta*wnl ) ; 

D1=[0;0;FC(1;1) /M(l;l) ) ; 

[yl] =lsim(Al, Bl, Cl, Dl, Force, t) ; 

%State  Space  Matrix  for  Equation  #2 
A2=  [0 , 1;  -  (wn2''2  )  ,  -  2*zeta*wn2]  ; 

B2=(0;FC(2:2) /M(5:5) ] ; 

C2=[l,  0;0,  l;-(wn2''2;  , -2*zeta*wi-2  ]  ; 

D2=[0;0;FC(2 :2) /M(5:5)]; 

[y2]=lsim(A2,B2,C2,D2, Force, t) ; 
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%State  Space  Matrix  for  Equation  #3 
A3=  [0,1;-  (vm3  ''2  )  ,  -2*zeta*wn3  ]  ; 
B3=[0;FC(3:3)  /M{9:9)]; 

C3= [ 1 , 0 ; 0 , 1 ; - (wn3  ^2 ) , -2*zeta*wn3 ] ; 
D3= [0;0;FC(3 : 3 ) /M ( 9 ; 9 ) ] ; 

[y3 ] = lsim(A3 , B3 , C3 , D3 , Force, t ) ; 

%Couple  the  Equations 


XIDIS 

= 

ull*yl( : ,1) 

+  ul2*y2 ( : , 1) 

+  ul3*y3 ( : , 1) ; 

XIVEL 

ull*yl( :  ,2) 

+  ul2*y2( :  ,2) 

+  Ul3*y3 ( : , 2) ; 

XlACC 

= 

(ull*yl( :  ,3) 

+  ul2*y2 ( :  ,  3  ) 

+  ul3*y3 ( : ,3) 

)  /g; 

X2DIS 

— 

u21*yl( :  ,1) 

+  u22*y2( :  ,1) 

+  u23*y3 ( :  ,  1)  ; 

X2VEL 

= 

u21*yl( :  ,2) 

+  u22*y2 ( : , 2 ) 

+  u23*y3 ( : , 2)  ; 

X2ACC 

= 

(u21*yl( : , 3) 

u22*y2  (  :  ,  3  ) 

+  u23  *y3 ( ; , 3 ) 

)  /g; 

X3DIS 

u31*yl( : , 1) 

+  u32*y2 ( ; , 1) 

+  u33  *y3 (  ;  , 1 )  ; 

X3VEL 

u31*yl( ; ,2) 

+  u32*y2 ( : , 2 ) 

+  U33*y3 (  : , 2)  ; 

X3ACC 

=: 

(u31*yl( : , 3) 

+  u32  *y2 ( : , 3 ) 

+  u33*y3 ( ; , 3) ) /g; 

plot ( t , XlACC ) , grid 

xlabeK'TIME  (sec )'), ylabel {• ACCELERATION  (g'’s)') 

title  (  'MODELED  THREE-DOF  FIXTURE  UPPER  TIER  ACCELERA'J. -iON '  ) 

gtext('5  Foot  Hammer  Drop  (200  g''s)’) 

gtext('Zeta  =  0.04') 

gtext( 'Upper  Tier  Weight  =  3500  lb') 

gtext( 'Lower  Tier  Weight  =  2700  lb') 

gtext( 'Anvil  Weight  =  4500  lb') 

meta  fig33.met 

pause 

%Upper  tier  velocity 
plot ( t , XIVEL ) , grid 

xlabeK'Time  (sec )'), ylabel ( 'Ve.ocity  -  ft/sec') 
title (' PREDICTED  3-DOF  MODEL  UPPER  TIER  VELOCITY') 
gtext( 'Wl=43 .5  lbs;  W2=57  lbs;  W3=109  lbs') 
gtext('Drop  Height=6  in;  Zeta=0.02') 
meta  threedof.met 
pause 

%Relative  Displacements 

plot(t, (X3DIS-X1DIS) *12, t, (X3DIS-X2DIS) *12) , grid 

title (' PREDICTED  RELATIVE  DISPLACEMENT  WITH  RESPECT  TO  ANVIL 

TABLE • ) 

xlabeK'Time  (sec )'), ylabel (' Displacement  (in)') 

gtext( 'Wl=43 .5  lbs;  W2=57  lbs;  W3=109  lbs') 

gtext('Drop  Height=6  in;  Zeta=0.02') 

gtext( 'Upper  Tier '), gtext (' Lower  Tier') 

meta  threedof.met 

pause 
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%Calculate  the  shock  spectra 
NF=120; 

SF=1; 

DF=1; 

Preq=[SF:DF:NF] ; 
vm=2*pi*Freq; 

F=[0;-1]  ; 

H=0; 

for  i=l :NF 

E=[0,l;-(wn(i)^2),0] ; 

GG=[-(wn(i)"2)  ,0]  ; 

[spec]=lsim{E,F,GG,H,XlACC,t) ; 
maxspec ( i ) =max ( abs ( spec ) ) ; 

end 

plot ( Freq, maxspec ), grid, xlabel {' FREQUENCY  (Hz) ') 
ylabeK 'ACCELERATION  (g''s)*) 

title  (  'MODELED  THREE  DOF  FIXTURE  ACCELERATION  SHOCK  SPECTRA’  ) 
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